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ABSTRACT 


The mechanics of alluvial channel bed-forns were investigated 
experimentally. Instrumentetion to measure flow velocity, pressure 
distribution and geometry of bed-forms was adspted or devised, end used 
4n obtaining measurements above and within actively translating, light 
waight sediment bedeforms in a laboratory flume. Velocities were meas 
ured with the hydrogen bubble technique, Pressures were sensed at 
twelve different positions within 9 bed-form, A series of smell diam- 
eter stainless steal wires, with a painted sesle and installed verti- 
cally in the flume were used te measure the geometry and te trace the 
bedeforms during their passage through the instrumented area. 

Preliminary investigations showed thet distortion of bedforms 
into three-dimensional configuration was caused ty side wail effects 
which retarded the flow and transport of sediments, A manual technique 
was devised to prevent bed-fora distortion end ail observations were 
made with tyoedimensional bed-form in the dune phase, Particle move- 
ment was in the pure bed-load mode. No suspended materials were 

The stability of all bed-forns was analyzed statistically and 
only those bed-forms which remained constant in length, height, and 
depth of submergence were analyzed. 

Flow shear was estimated from velocity profiles using rigid 
boundary techniques. 

The results showed that the maximum pressure gradient on bed-forms 
was maximum et locations of maximum bed-form surface slope and, hence, 


at locations of maximum particle entrainment, The minimum pressure 
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gradient occurred slightly upstream of the crest and corresponded with 
the location of zero particle entrainment. From the reattachment point 
the boundary shear increased in the downstream direction reaching & 
maximum at approximately 1/4 of the bed-form wavelength upstream of the 
crest and decreased to zero at the crest. 

The relative contribution by mean boundary shear and mean pres- 
sure gradient to particle entrainment showed that the mean pressure 
gradients exert finite forces upon the stationary particles on the sur- 
face of bed-forms, The effects of dilation upon particle concentration 
and friction in the surface layer were identified as critical in the 
evaluation of relative contribution to entrainment by different hydro- 
dynamic forces. 

Velocity profiles over a section of upstream face of bed-forms 
were found to be logarithmic, High velocity flow cell and a zone of 
turbulence were found to extend to the free surface of the flow above 


the trough region of a bed-form. 
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LIST OF SYMBOLS 


distance from hydrogen bubble wire to a point in the row of 
hydrogen bubbles 

regression coefficient 

constant 

volume concentration of particles 

constant 

depth of flow 

mean depth of flow above a bed-form 

mean depth of flow above crest 1 

mean depth of flow above crest 2 

diameter 
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is finer by weight 

minimum force required to disturb particie equilibrium 
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to disturb particle equilibrium 
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required to disturb the equilibrium of the surface layer of particles 
gravitational acceleration 

volumetric rate of bed-material entrainment per unit width of 
bed-form from segment 0.1L long 

volumetric rate of sediment flux per unit width of flume at the 
downstream limit of a segment 


volumetric rate of sediment flux per unit width of flume at the 


erest of a bed-form 
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vertical coordinate of bed-form surface at segment limit 
instantaneous vertical coordinate of the origin of the moving 
reference frame 

instantaneous vertical coordinate of a point on the bed-form 
surface in the stationary reference frame 

local hydraulic pressure gradient 
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in the moving reference frame 
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profile 


mean height of bed-form 1 

mean height of bed-form 2 

mean reference height of a bed-form 

vertical coordinate of a point on the bed-form surface in the 
moving reference frame 

wavelength of bed-form 

mean wavelength of bed-form 
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mean observed pressure 

mean reference pressure at the origin of the moving reference frame 
dimensionless pressure 
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time 

temperature 
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friction velocity 

mean instantaneous flow velocity 
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velocity,. all-profiles 

mean velocity based on flow discharge and average depth of flow 
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mean velocity of bed-form translation 

mean velocity of flow above the crest of upstream bed-form 
total volume occupied by solid particles and fluid 

body force due to gravity on all particles contained in a volume 
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body force on a particle due to gravity 

instantaneous horizontal coordinate of the origin of the moving 
reference frame 
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CHAPTER IT 
INTRODUCTION 


1.1 Introductory Remarks 

A channel is said to be alluvial if at some stare of flow its 
non-cohesive boundary materials are entrained, transported, and deposited 
by the flowing water. The hydraulics of such channels are very compli- 


cated because the resulting water-sediment sompisx constitutes a two 


flow, establishes the geometry of the channel boundaries, and moves the 


sediments by different modes of transport. The ability to predict these 


naturally developed features Bae ts alluvial system often determines the 
success or failure of a water resources project in terms cf economical, 
environmental, and even sociological PS reise eetaene™ 

For nearly one hundred years engineers and geologists have sought 
_ to understand the mechenisms of flow in alluvial channels, yet relatively 
little has been achieved in the analytical formulation of the phenomena, 
The progress toward a comprehensive solution of the alluvial flow problem, 
measured by the stage of development within the context of fluid mechanics 
research has been reviewed by Engel (1969). Using a "knowledge develop- 
ment" progress scale proposed by von Karman in 1934 Engel suggested that 
the scientific inquiry into this problem has progressed from the initial 


stage of empirical laws into the intermediate stage of knowledze develop- 


ment. Many variables controlling the phenomena have now been identified 


1 
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by systenatic analysis of the empirical information. Some qualitative 
and quantitative aspects of the interdependence of various variables have 
been identified and expressed in functional relationships among different 
sets of numerics (dimensionless parameters). Some of these relationships 
may now be used to predict practical flow behaviour within tolerable 
accuracy limits. In addition entrainment of particles, their transpor- 
tation along the bed or in suspension, channel boundary stability (bed- 
wave and meancer formation) and channel roughness have been identified 
and recognized as integral components of the overall problem. The lim. 
ited progress toward the ultimate goal of establishing a comprehensive 
alluvial channel flow theory can be partially attributed to the failure 
on the part of researchers to explore and understand the basic physical 
features of the phenomenon. In the apsence of this information many of 
the proposed theories have been based on assumptions which most often 
were incompatible with the natural state of flow. Hence such theories 
can only be regarded as approximate representations of the natural laws. 
It appears, therefore, that the logical research direction toward 
the establishment of reliable theories for solution of problems of flow 
in alluvial channels should commence with a systematic experimental 
exploration of the physical mechanisms responsible for the entrainment, 


transport, and deposition of solid particles by the flow. 


1.2 Bed-Waves 

Perhaps the least understood aspect of alluvial flow phenomena 
concerns the generation and occurrence of trains of periodic bed-forms 
on the bed of alluvial channels which form from non-cohesive bed materials. 
The importance of bed-forms from a practical point of view has long been 


recognized because their presence or absence to a large degree determines 
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the roughness and hence the flow resistance in the channel, Furthermore, 
the bed~forms are related to the transportation of material along the bed. 
Two general directions of approach, distinguished by the purpose 
for which the information was sought were employed in the study of bed. 
forms, Although concerned primerily with engineering design information, 
the empirical studies recognized at least six distinct bed.form confi- 
gurations each associated with different conditions (regimes) of flow in 
the channel. This suggested that different types of bedforms were genere 
ated and maintained by distinct dynamic flow systems and raised the possi- 


bility precluding the existence of an universal theory for #11 bed-forms 


(Blench, 1969). 


stochastic techniques or analytical models and aimed toward improvement of 
understanding of physical laws sovarning the bed-form processes have been 
far less successful in Hdadhone their objectives, From the results of 
these complimentary studies, however, two conclusions were drawn and were 
used to evolve the present investigation: (1) improvement of understanding 
of bed-wave mechanics can be expscted only from detailed analysis of kins- 
matic and dynamic interaction between the fluid and bed-material in 
alluvial channels (Kennedy, 1969), and (2) systematic verification of 


physical facts by experiment should procede analytical formulation of 


theories (Rivlin, &s reported by Davidson et al. 1969). 


1.3 Purpose, Objectives, and Scope 
fhe purpose of this study was to investigate experimentally the 


A4nteraction between the flow and the non-cohesive bed-material particles 
in alluvial channels. It was intended to extend and improve the under- 


standing of the physical. mechanisms responsible for the orderly scour- 
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transportation-deposition of sedimentary materials in the presence of 
translating bed-forms, 

The investigation was entirely experimental, Unlike in previousely 
reported studies, actively moving alluvial bed-forms were employed for 
observétion of various flow and boundary geometry parameters, 

The main objective of this investigation was to verify the general 
hypothesis that entrainmert of stationary bed particles from the surface 
of active bed-forms in alluvial channels was the result of a combination 
of hydrodynamic forces which are inherent in the equations of fluid 
motion, The specific objectives were (1) to develop reliable experimental 
instrunentation and methodology for studies of moving alluvial bed-forms 


in 4a lsaherstary flume, 
Ae _— ete ST oO ee =e ~ 


y ume, (2) to ebtain experimental data and to establish 


the characteristics of various flow, sediment response, and bed-form 
geometry parameters, and (3) to evaluate and correlate the relative 
magnitudes of the spatial distribution of available forces for the entrain- 
ment of bed particles from the surface of bed-forms. Special attention 
was focused upon the effects of pressure variability upon the sediment 
perticle entrainment by the flow. 

A number of limitations were imposed upon the scope of investiga- 
tion. Since the aims of the present investigation were of fundamental 
rather than practical significance, to facilitate the experimental inquiry 
and to simplify anelysis, idealized (light weight, spherical shape ) 
bedenaterial was used and only two-dimensional bed-forms in the tranquil 
flow regime were studied, Flow velocity measurement technique restricted 
transport of the sediment to the pure bed-load mode. 

No attempt was made to utilize the experimental results for 


improvement of mathematical models of bed-form mechanisms, The investiga- 
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tion was intended to demonstrate physical facts about the finid-sediment 
interaction phenomena and to establish a precedent for a new type of 


experimental approach as a guide for future studies, 


Suensivx— tra of? J0e0R sjoe? fon tepety: atertonnentb ad fohownt eau a 


ers. 





1a eae won 4 161 Saehesets « de bigaiey oF tut enemeeesiy 


tarhvte erptat-agl autem § ce woeemgee ies 





= * 
.* .* * 
ae td | . 
» ‘ 
: 7% ” ¥ + . c 
v «- = - « . - ‘ ° > - ¢ 
“ - $ - - ~ ~ 
i * 
€ * » a F S % % ° a ee - * 
* 
= - * * S - 
% ’ . « 4 a tr “Ae 
* 
« é > + + 
. e * a 
7 
. 
-~ ~ “ a* * ‘ « -< . 
: ~ ais 
7 
é = © 4 i sd . 4 "4 , - 
’ = © : .¥ i » 4 =_4 a 
~ * “ a * 
. : Sig gia Re a0 9 +: hemi i 





CHAPTER IT 
LITERATURE REVIEW 


2.1 General 

The interest in sedimentary bed-forms can be traced to the last 
century when Darwin (1883) published the results of experiments of 
ripple formation in oscillatory flow. Since that time numerous other 
researchers approached and investigated the problem from different 
directions, Attempts have been made to explain the bed~form phenomena 
through classical mechanies concepts, water-sediment interface instabil- 
ity considerations, kinematic wave models, by dimensional analysis con- 
siderations and through statistical techniques. And, although consid- 
erable insight has been ebteined from these studies, none of the above 
approaches has been entirely successful. in explaining the bed-form 
phenomena. 

According eater cer (1971) three separate but interdependent 
factors must be considered in the analysis of the bed-form problem: 

(1) The flow of water above the bed-form. 

(2) The geonetry of the bed-form. 

(3) The transportation of bed particles over the bed-form. 
A bed-form model, therefore, must incorporate three separate relation- 
ships: | 

(1) The transport of sediment and fluid flow over the bed-forn. 

(2) The fluid flow and the geometry of the bed-form profile. 


(3) The transport of sediment and the geometry of the bed-form 
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profile, 

Of the three relationships only the third, concerned with the 
continuity of sediment motion has been formulated exactly. However, in 
various analytical models, the other two relationships were based on 
various simplifying assumptions. As pointed out by Mercer (1971) in 
the development of the bed-load transport relationships a number of 
factors, such as pressure gradients, bed-slope and boundary layer de- 
velopments, which may become important in bed-forms have generally been 
ignored. Similarly the analytical models generally employed simplified 
models to represent the flow over the bed-form shapes, For example, 
Exener (1925) considered a simple one-dimensional flow above the bed 
and assumed the local flow velocity to be inversely proportional to the 
local depth of flow. Both Anderson (1953) and Kennedy (1963) assumed 
potential flow. 

Although the various eee cad models are found to be highly 
instructive, they nevertheless have failed to explain satisfactorily 
the mechanics of bed-forms. The main reason for this failure can be 
attributed to the fact that the assumptions used in developing the re- 
lationships between the flow and the bed-form geometry, and between the 
flow and the transport of bed-load used in the models did not fully 
consider the dynamic and kinematic aspects of interaction between the 
fluid and sediments and were inconsistent with the actual conditions 
eccurring in the two-phase systems. 

A review of previously reported experimental investigations 


related to bed~form mechanics is presented in the subsequent sections. 


2.2 Experimental Approaches 
The majority of all experimental investigations concerned with 
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bed-forms and reported in literature have been oriented toward explain- 
ing the gross features of the phenomena, especially as it relates to the 
roughness characteristics of the channel. Only a very small number of 
investigators have concerned themselves with the flow and sediment 
response conditions within the macro-environment of individual bed-forms. 

Furthermore, literature survey disclosed that, because the prob- 
lem is very difficult experimentally, no attempts have ever been made to 
study the phenomena using active bed-forms, This has been recorpnized 
by Vanoni and Hwang (1967) who observed that experimental measurements 
are difficult to make under natural conditions when bed-forms translate 
and the presence of instruments may modify the sediment surface and 
bed-form patterns. 

Because of these difficulties all previously reported investi- 
gations were carried out using stabilized bed-form shapes. The stable 
bed-form approach, however, only approximately represents the true con- 
ditions of flow above the bed-forms. Whether the bed-form shapes are 
constructed from solid materials or by spraying with adhesives, the 
natural bulk characteristics of the granular materials are no longer 
properly represented, In the absence of particle movement the observed 
velocity profiles, likewise, may not represent the true flow velocity 
conditions which exist in the presence with sediment transport. Never- 
theless, these investigations have provided considerable insight into 
the Gow characteristics which exist above individual bed-forms and 
offer the only experimental background for a highly complex problem in 
hydraulics, 
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2.3 Hydrodynamic Force Systems 

2.3.1 Introductory Remarks 

The occurrence of bed-forms in alluvial channels is dependent 
upon an orderly pattern of scour, transport, and deposition of sediment 
particles, Under equilibrium flow conditions, the flow induced forces 
which act upon the sediment and the response of the sediment particles 
to these forces are mutually compatible in that relative to the geometry 
of the bed-form a steady state of flow, sediment entrainment, transport, 
and deposition is observed in the mean above the bed-form. 

Previous experimental results and general observations of 
particle entrainment, transport, and deposition above bed-forms have 
indicated that the forces responsible for the sediment responses vary 
spatially above the bed-form. Furthermore, there is evidence which 
suggests that these fluid forces may derive from different origins in- 
herent in the equations of motion. Hence, boundary shear, pressure 
gradient induced seepage, and turbulence, as well as particle impact 
and gravity may all contribute to the detachment and transport of 
sediments. 

2.3.2 Boundary Shear Stresses 

Traditionally, boundary shear stresses have been considered as 
the main cause of particle entrainment. This, perhaps may be traced 
back to the work by Shields (1936), who successfully established incip- 
jent motion relation for flow over flat bed with shear as the main 
parameter. The Shields! criteria, however, although fundamentally sound 
in the particular context in which they were developed, do not represent 
the conditions existing above natural bed-forms, 


Laursen et al. (1962) investigated flow of air over schematic 
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triangular shaped bed-forms placed symmetrically at top and bottom of 
a closed conduit. eee shear was evaluated by the Preston technique. 
The measurements showed that shear increased almost linearly downstream 
from zero at reattachment point and was maximum at the crest of the 
Deere tra 

Raudkivi (1963), in a classic papsr to-date, measured velocity 
profiles and evaluated bound&éry shear due to flow over a bed-form 
shaped from galvanized sheet metal, His main conclusion was that entrain- 
ment and transport of grains were not unique functions of the temporal 
mean shear or mean velocity, The shear increased from zero at the 


reattachment point and was maximum at the crest. The critical shear stress 


be half-way up the upstream slope of the bed-form, which is considerably 
dewnstream from the position where actual entrainment of particles tock 
place, 

Khanna (1970) also investigated boundary shear of air flow over 
soft wood bed-form shapes. However, no flow separation occurred at the 
erest and shear was found to be finite in the downstream direction every- 
where above the bed-form. 

2.3.3 Pressure Distribution Effects 

Pressure and pressure gradients have long been suspect as factors 
in the erosion and transportation phenomena, Blench (1969) has reported 
observations of piping failure of antidune slopes. Simons (1962) noted 
increase in total transport of sand in alluvial channels due to surface 
water waves. Posey (1963) suggested that when high and low pressure 
gsones exist in the vicinity of bridge piers in porous bed-material the 


resulting seepage enhances the removal of the surface material. 
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Pressure distributions over the stabilized triangular bed-forms 
have been measured by Laursen et al. (1962) who found that a systematic 
difference existed between the pressure difference for different 
surface roughness, 

The pressure measurements by Raudkivi (1963) established that the 
maximum pressure occurred at the eine of the separation streamline re- 
attachment on the upstream face of the bed-form. This was confirmed by 
Vanoni and Hwang (1967), and Rifai and Smith (1971). 

The position of the minimum pressure intensity on the surface of 
the bed-forms has not been as well established. On the triangular bed- 
forms reported by Laursen et al. (1962) the minimum pressure appears to 
be at the crests of the bed-forms. The minimum pressure from Raudkivi's 
-(1963) measurements occurred slightly upstream of the crests; however, 
Vanoni and Hwang (1967) and Khanna (1970) found that the minimum pressure 
was on the lee face of the bed-forms. 

Generally, the pressure distribution curves established by 
measurements over the stabilized bed-forms downstream from the reattach- 
ment point were found to be concave downward to approximately half-way 
between the reattachment and crest positions. The position of pressure 
distribution curve inflection point, however, has not been specifically 
identified by any of the authors. 

The contribution of pressure gradients established by the flow 
above granular material bed-forms have not been explored to-date. 
However, a theory was proposed by Luque (1967) which takes into account 
the seepage flow through the granular media and the forces which act on 
the grains near the bed surface, The guiding principle of this theory 
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“the rate of change of momentum in ie turbulent main flow 
with the transport of momentum across a turbulent boundary 
layer towards the bed surface and with the further transmission 
of stress through the solid and fluid components of the bed." 

This theory recognizes that pressure gradient induced seepage 
may exert considerable drag forces upon the bed surface particles, It 
indicates that a bed of loose granular sediment is severely eroded at 
locations where the momentum of the surface flow strongly changes, and 
attributes ripple occurrence in sand bed channels to strong pressure 
gradients in the vicinity of stagnation points of standing vortices 
formed behind discontinuities of the bed surface, 

| It must be noted that this theory was questioned (Davidson et al. 
1969) on the grounds that previously neither Rouse (1940) nor Laursen 
(1952) found any evidence to indicate that the seepage forces were 
| Significant. Yet, a review of these papers disclosed that pressures 
were not considered by these two investigators. 

The effects of seepage force on bed particles comprising the 
top layer of a stream bed have been investigated ane. Pimentel 17 tby 
Martin (1971), These experiments were concerned with seepage normal to 
the bed and disclosed that seepage force on the top grains was approxi- 
mately one-half that on grains several layers within the bed. 

An extensive review of the effects of seepage into and out of 
the granular bed has been presented by Watters and Rao (1971). Their 
main conclusion was that seepage modifies the sublayer and velocity 
profile characteristics near the sediment bed. 

Effluent seepage from the bed was found to decrease the drag 
upon the particles regardless of their position on the bed. It increased 


the lift upon a plane bed particle but decreased the lift if the particle 
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was above the general bed level. | 

The opposite result was found to hold for influent seepage into 
the bed. 

2.3-4 Turbulence Effects 

Turbulence has long been recognized as a potential contributor 
to the causes of sediment erect and transport. The qualitative 
aspects of turbulent flow upon the detachment of particles have been in- 
vestigated by Sutherland (1967) who established entrainment hypothesis 
based upon the concept of turbulent eddies disrupting the viscous sub- 
layer and impinging directly onto the grain surface. The effects of 
dynamic tubulent pressure fluctuations upon the bed particles have been 
recognized also by Kontratev et al. (1962). 

Turbulence intensity measurements over stabilized bed-forms have 
been made first by Raudkivi (1966). The results have shown that maximum 
turbulence energy is reached in the vicinity of the reattachment point and 
decreases rapidly in magnitude toward the crest. The results presented 
by Khanna (1970), however, indicate that the highest longitudinal inten- 
sity of turbulence occurred in the trough region of the bed-form. These 
results, however, may have been influenced by the fact that flow did not 
separate at the crest, 

Rifai and Smith (1971) observed turbulence characteristics and 
found that the macroscale of turbulence is related to the height of the 
bed-form elements. 

More recently, Williams and Kemp (1971) established that close 
to flat beds the longitudinal velocity Pie erinurtons had a high positive 
skewness caused by a number of instantaneous very high velocities. 


However, to-date no technique for evaluating the forces upon 
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sediment particles due to turbulence fluctuations above bedforms has 


been evolved for pratical. applications, 
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CHAPTER III 
LABORATORY EQUIPMENT AND INSTRUMENT ATION 


3el Introduction 
Because of the uniqueness of the experimental environment and 

the absence of & precedent for guidance in selecting apparatus for 
usssuronent of flow and bed geometry parameters, part of this investiga. 
ticn was concerned with the development and evaluation of instruments and 
technigues of measurement. Conplete description of laboratory equipnsnt, 
bedenmsisrial, and instrumentation systems are presented and discussed in 
this chapter. 

| The discussion ealatonato instrumentation falls into three broad 
categories, First, the environnental characteristics relevant to the 
measurement of the various paremeters are discussed, Alternative cone 
cepts are examined under the restrictions imposed by the test conditions. 
Next, under soparate headings for each of the measured parameters, the 
respective probe~data acquizition systems are described, Finally, system 
evaluation study results, to show expected measurement accuracy end 


compatibility of each system with expsrimental environment are presented. 


3.2 Laboratory Equipment 

32.1 Flume 

The expsrimental investigation was carried out in a rectangular 
flums which was located on a special platform 7 feet above the main 


laboratory floor in the Hydraulics Laboratory of the University of Alberta, 
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The flums was 119 feet long, 16 inches wide, and 16 inches 
deep, The sides and the bottom of the flume were of 3/8 inch thick 
plexiglass panels whick were supported by continuous longitudinal two 
inch aluminum en¢les, The vertical joints between adjacent panels of the 
plexiglass walls were at 8 fect centers and were filled with carefully 
smoothed silicone sealant to eliminate swrface roughnesses and make the 
flune watertight. 

The headbox at the upstream end ef the flums was 2&8 inches wide 
and 19 inches deep. A six inch radius transition at the sides and the 
botten guided the flow from the headbox inte the flums. At 3 feet from 
the downstream end, the flume was stepped down to a depth of 15 inches 
and a gradual transition was previded into a 4 inch diameter pipe. 

The fluns wes supported cbowa the floor by H-frames made of one 
inch diameter pipes, Laterel end lercitudinal crossebracing was provided 
by emell diameter galvanized steel cables tensioned with turnbuckles. 

3o20% Flow System 

The flow system eanloyed in this investigation was of the ree 
cirveuleting types in which the discharged wateresediment mixture was 
pumped back into the headbozx of the flume, A schematic of this system 
is shown in Figure 3.1. 

Addittenal equipment used in the system included the following: 

(1) Vertical-centrifugal, five inch diameter pump of the closed 
4npoller type, driven by ® constant (1150 rpm) speed, 7.5 horsepower 
electric motor. 

(2) Four inch gate valve for the control of the flow. 

(3) An eloctronagnetic flow moter (rated at 0.01519 millivolts 


per U.S, gallon por minute) with a recorder for volumetric discharge 
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measurement te en accuracy of 0.01 cubic fest per second. 

(4) Four inch polythene pipe to convey the flow-sediment mixturs 
to the headbox of the flume. 

(5) Plywood plate 24 inches long, 16 inches wide, to suppress 


surface waves at the upstream end of the flume. 
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Because this investigation was concerned with fundamental rather 
than practical. aspects of the alluvial flew problem, idealized bed. 
material wes used in the experiments, This simplified the interpretation 


of results, exaggerated some of the processes, and facilitated some of 


White spherical polystyrene particles were chosen for the experi- 


e 


rents, A nacresphotograph ef the particles is given in Figure 3.2, 


Ss “Fan 
AN 





Figure 3.2. Photograph of Sediment Particles 
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The photograph was taken at the conclusion of the testing program, and 
the two non-spherical particles seen left of Panter on the photograph 
indicate that some grinding of the particles occurred during the movement 
of the sediment through the pump. | 

|—--3.362 Specific Gravity 

The average specific gravity of particles was found to be 1.041. 
It was determined by general procedures of the American Society for 
Testing of Materials standard ASTM D85%4.-58 (1965) for soil particle speci- 
fic gravity. However, it was necessary to reduce the specified particle 
drying temperature to 30 degrees Celsius because temperatures above 40 
degrees Celsius caused discolouration of the polystyrene. 

The results of tests to determine variations in the specific gra- 
vity of individual particles are summarized by the histogram shown in 
Figure 3.3. The procedures used to determine the distribution of specific 
gravities of the sediment and the results of four individual tests are 
given in Appendix A, 

3.33 Grain Sise Analysis 

The grain size distribution of the sediment was determined in 
accordance with the American Society for Testing of Materials standard 
ASTM D422-63, The results are summarized in Figure 3.4. Approximately 
97% of the particles range in size from one to two millimetres. The 
median particle diameter Deo (the size of sediment of which 50% is finer 
and 50% is larger by weight) is 1.4 millimetres. 


The gradation coefficient 


4s found to be 1.26. The uniformity coefficient 
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Figure 3.3. Histogram of Specific Gravity Distribution 
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was calculated to be 1.27. According to Jumikis (1962), a uniformity 
cosfficient of less than five represents a very uniform grain size dis- 
ceUWtion. 

3.3.4 Permeability and Porosity 

The permeability characteristics of the sediment material used 
in this study were determined using the equipment and procedures 
described in the American Society for Testing of Materials standard 
ASTM 2435-6557 for permeability of granuler soils under constant head. 
The average coefficient of permeability (at 20°C), determined from five 
separately prepared samples, was 0.121 centimetres per second. 
Sinultansously determinad sample average porosity and volumes concentration 
values were 6.402 and 0,598 respectively. 


3.3.5 Angles of Initial Yiold and Repose 





A plane surface underlain by loose granular particles can be 
built up in the gravity field only up to 4 certain critical slope angle, 
c@lled the ati of initial yield ¢4. When this angle is exceeded, 
avalanching takes place and the slope assumes & smaller angle dy 
commonly called the angle of repose (Allen, 19702). 

The angles $d; and >, of the polystyrene particles in water used 
in this investigation wore determined experimentally. A closed glass 
container, approximately one litre in velumo, with opposite side walls 
plane and psraliel to each other, was filled half-full with the granular 
perticles, Water wes added to the brim of the container which then was 
closed excluding all air from within. The container was turned on its 


side, and moved slightly to and fro to establish a flat interface between 
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the surface of particles and the water above, It was then tilted slowly 
from one end in one degree increments until a general renanti ot parti- 
cles was observed, The angles of tilt of the original surface to the 
horizontal at the initial movement and the final accumulation of parti- 
cles, ds and >. respectively were recorded, From ten observations, the 
averages were found to be by = 31° and >. = 26°, 

It was found, however, that $4; was affected by the method of 
the surface. particle deposition. If a surface of particles deposited at 


g, was tilted in the opposite direction, 4 was found to be 29 degrees. 


324 Experimental Instrumentation 
3.4.1 General Considerations 


Three distinct aspects of the experimental environment pertinent 
- to the development of instrumentation systems for measurement of flow 
and sediment bed parameters were identified: 

(1) Time dependent (unsteady) nature of the two-phase flow 
above bed-forms,. 

(2) Sensitivity of the natural equilibrium between the flow 
and the granular bed-material,. 

(3) Low absolute magnitudes of velocities and pressures near 
the sediment bed. 

Under established equilibrium conditions a bed-wave retains a 
constant profile and a steady but varied mean structure of flow above. 
However, because the bed-wave and its Remctatad flow structure trans- 
late, the entire water-sediment system becomes unsteady when viewed from 
a stationary frame of reference. 

To determine experimentally the mean steady state flow characte- 


ristics associated with a translating equilibrium bed-wave, measurements 
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of the various flow parameters at a number of different locations in the 
flow system are required. To eliminate the unsteadiness factor, which 
affects the positioning of the sensor probe in relation to the bed-wave, 
two alternative approaches were considered, In principle, a system of 
sensor probes installed on an auxilliary vehicle traveling at the speed 
of the bed-wave would have served ine purpose. However, this concept 
was found to be impractical because measurements were to be made within 
the mass of bed-wave particles, When a bed-form translates, only a 
small number of grains on its surface is in actual motion at any time. 
Movement of sensor probes through the stationary mass of granular 
particles would have disturbed the stability of boundary materials and 
consequently the equilibrium of the entire system. 

To eliminate the unsteadiness effect all sensor probes were 
"permanently installed in the flume. The instantaneous positions of the 
probes relative to the moving bed-wave were determined by observing the 
instantaneous horizontal positions of the bed-ewave crest, 

The major disadvantage of this method was that observation of a 
flow parameter could not be repeated at the same, relative to the 
bed-wave, position with the same sensor probe. On the other hand, a 
sensor probe with an analog output could provide an uninterrupted 
picture of the spatial changes of a measured parameter as long as its 
relative position in the moving reference system was Sos 
monitored. 

It has long been recognized in instrumentation engineering that 
the introduction of a finite size sensor into a flow modifies the field 
of the parameter the sensor is to measure, As a result of this pheno- 
mena two consequences of direct interest in this investigation were 


identifio?t; (1) the probe output may not represent the true value of 
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the measured parameter, and (2) the presence of the probe may destroy 
the natural equilibrium between the flow and the sedimént boundary, 
The new equilibrium, compatible with the combined flow-sediment-probe 
interaction system may be totally different from the natural system 
under investigation. 

In most fluid flow measurements, the flow field modification 
effects upon the sensor may, for practical purposes, be eliminated by 
calibration of the measurement device irae model environment known to 
represent the conditions in the prototype. However, because this study 
was concerned with measurements in an unknown environment, the modelling 
approach could not be employed. } 

In principle, at least three alternative approaches may be used 
to eliminate the flow field modification effects due to sensors on their 
output and on the boundary stability: (1) reduction of the physical size 
of the probe elements to zero in the limit, (2) matching of the physical 
properties of the probe elements to those of the medium in which the 
measurements are to be made, and (3) employing remote sensing techniques. 

None of these alternatives is technically feasible at present. 
Remote sensing techniques must be perfected by further developments, 
matching ot pryettal properties of sensor probes to those of the medium 
is impossible even in single phase flow environments, and sensor probes 
have not been reduced to zero size. Therefore, the compromise solution 
was to use sensor probes with minimum physical dimensions, 

Simple preliminary tests were carried out to ascertain the 
maximum permissible size of measurement probes that would not disturb 
the natural water-sediment boundary equilibrium, In the anticipated 


range of average flow velocities from approximately 2.5 to 4.5 inches 
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per second, cylindrical shapes protruding vertically through the 
sediment boundary into the flow did not disturb the appearance of the 
sediment bed if the Reynolds number, based on the diameter of the 
cylinders and average flow velocities did not exceed approximately 350. 

The need to measure small absolute values of velocity and pres- 
sure may not have been unique in these experiments, Nevertheless this 
requirement precluded the use of standard laboratory equipment (Pitot 
tubes and prpeller meters) for the measurement of velocities, Simi- 
larly it restricted the choice of pressure measurement equipment to the 
most sensitive, low range commercially available units, 

3.4.2 Bed-Wave Profile Measurement 

To measure the geometry of the bed-waves, a series of nine 0.012 
inch diameter stainless steel wires were employed, Alternate half-inch 
yer along the length of each wire were painted red to provide a re. 
ference scale, The wires were installed vertically between the bottom 
of the flume and a support bar above, They were spaced at two or four 
inch intervals along the longitudinal centre line of the flume. Details 
of a depth-wire and its support are shown in Figure 3.5 (a). At the 
bottom the wires were held in position by a threaded brass plug screwed 
flush with the top of the plate, The wire passed through a hole at the 
centre of the plug. The wire was prevented from pulling out from the 
plug by a knot tied on its end. | 

At the top, each wire was attached to a small tension spring 
which was held in position by a hook installed on the underside of the 
supporting bar. 

Because the scale divisions on the depth-wires were one half 
inch the accuracy of bed-wave profile measurement depended upon the 


observer's ability to interpolate the position of the bed level within 
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the one half inch division. The precision of the depth-wire technique 
was ascertained by observations of bed levels on residual bed-forms re- 
7 maining in the measurement area after the generating flow was stopped. 
All other conditions such as lighting, observer's position, were main- 
rained identical as in the Subsequent experiments where the bed-form 
geometry was measured, 

304.3 Velocity Measurement 

In view of the considerations in Section 3.4.1 a relatively new 
method of velocity measurement--the hydrogen bubble technique--was 
selected for this investigation. 

Detailed expositions of the hydrogen bubble techniques have been 
given by Clutter (1961), Schraub (1965), and Davis and Fox (1967). In 
the discussion that follows only those aspects of the technique which 
are pertinent to the present study in the sediment transport research 
area are considered, 

The hydrogen bubble method is primarily a visual method which 
utilizes the electrolysis of water at a very fine wire to introduce 
hydrogen gas into the flow. The hydrogen bubbles produced at the wire 
are swept away by the flowing water and act as tracers which can be made 
visible by lighting and hence can be photographed. If the electrodes are 
subjected to pulsating applications of electrical energy, rows of hydro- 
gen bubbles are produced. Since te frequency of pulsations can be 
controlled and measured, and the distances between adjacent rows of 
bubbles can be scaled from photographs the velocity at a point in the 
flow can be calculated by multiplying the true separation distance be- 
tween two adjacent rows of bubbles by the frequency of pulsation. 


The hydrogen bubble method possesses a number of advantages over 


vs PA ae | ns 


cat an 
avpimuiseds evtw-2qgel) oft to cobetomtg ae? aokeleth donk Lied ewe edt 
- Qeri@lebec Leublest ao tlevel bed Yo gyabtevsseds ve bon tataenss ter 
.netyota 2tw wor antiacenes #42 Tar) s afar toneewtsee ofa nf yittortonge 
-tfan ecaw .Woattiz e'vavoedo ,~ottaell ef dete @nolstbuce seddo TIA 
wtelebed.ai? ai10/ Ina iteqxe Jrecpecdsa eit mf aa Laokineh! baakay 
sbeissea saw rideneny 

inawer ure fener €. ¥ 

s I,u.f aolvose of sonktacehdensos af? Yo welv al 
’ inet afdcid vepoxbyr! siij«e4enengenem phloolévy to baties 


waivegivewtnt sit? «02 betoeles 


teed evad gseupiadise?s aficd neyotDiy! ert to eee lirsogme be Liaeted 
nl ,(§0@i) xo% dae etvad fe (ede!) deweded (ECL) retanid yw geeks 
é * exotd elooe awuifcl gad? nofeensetS ed? 


doyseteat s1ragtnes) Jeeniooe ate gt 4 ute tener ei? of Sneaniiseg ers 
| »bereb iene wra So Th 

iw -bopren Layety 6 vii ree boijen @elddud sesorigt eT 
scrmoysnl of exiw eat? ts07 & 74 suitiad 20 siaytoroele ed? eealhl itp 
etiw of) de heovbow aeldcud coporhyd eft - yell sdf oft? esy negotiyd 


ahem @0'aas-reoliw tioords ca tne bis tetew aeivell ad? wa ie 'qews ous 


ora tebortants ait LI .betwersoteig ed wae ener dae antdiigi ww aldiaky 
<othyd Yo awor .yuaernn LaatesouSa ‘to eattiaetigqgs aR ot 
—— eee. bed te le eg Pisce Pegg 


— 





28 


Depth Wire Hydrogen Bubble Wire 


: Teflon wedge 








Ficure 3.56 
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Sediment Depth Wire and Hydrogen Bubble Wire Installation 
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other conventional metheds of velocity measurement. First, the electrode 
wire used to preduce hydrogen has a very small dienmeter which has little 
effect upon the velocity field in its vicinity. The maximum Reynolds 
number of the tungsten wire used in these experiments was approximately 5. 
This was considerebly below the Reynolds number valucs at which the wire 
wake effects upon the sediment bed become noticeable, Secondly, hydrogen 
bubble photographs provide information for a simultaneous and continuous 
measurement of velocity over the entire depth of flow--a condition that 
cannot be realized by conventional point velocity measurement methods. 
Third, it provides visual information about the flow structure. Fourth, 
it can be used to measurs very low velocities, with the lower limit de. 
termined primarily by the bouyancy effects of the hydrogen bubbles in 
water. 

A major disadvantage of the method is that it yields instantanecus 
rathsy than mean velocities of flow. Statistical means can be employed to 
obtain mean velocity values if sufficient number of instantaneous records 
are taken, however, the reduction of date is extremely time consuming. 
Furthermore, transverse movenent of water from the plane of photograph 
cannot be detected, 

The apparatus to produce and observe hydrogen bubbles in water 
consisted of the following: (1) electrodes, (2) electrical apparatus for 
energizing of the electrodes, (3) illumination, and (4) photographic 
equipzent, 

Ths neg@tive electrodes, at which hydrogen bubbles were generated, 
wore three sepiratse 0,002 inch diameter tungsten wires, They were in- 
stalled perpendiculer to the direction of flow in the vertical plane of 
the longitudinal centreline of the flume, and were spaced at four inch 
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Esch wire was held tautly in position by an arrsngement shown in 
Figure 3.5 (b). A slightly tapered hole of small diameter was drilled 
through the centre of a 1/4 inch diameter threaded plug. An end of a 
tungsten wire was passed through the hole and jamzed against its side by 
a watching teflon wedge. The assembly then was screwed inte a threaded 
hole in the bettom plate. 

At the top, the wire passed through &® slightly tapered 1/ inch 
diameter hole in a 2% inch wide plexiglass bar, It wes tensioned by hend 
and held against the side of the hole by a matching conical teflon wedge, 

A 3/32 4nch diemeter aluminum welding rod was used as the positive 
electrode, It was supported from a plexiglass bar, spanning the sides of 
the flume, and was located approximately six inches and on the far side 
from the tungsten wires from the observation window. Depending upon the 
depth of flow the total submerged length of the anode ranged from 7 to 9 
dnches, To minimize tho disturbance of the flow by the rod, its submerged 
end was bent into an L-shape and the end leg was positioned parallel te 
the flow. 

The electrical system to produce hydrogen bubbles is shown in a 
schematic diagram in Figure 3.6, It consists of two groups of components, 
each serving a different function: (1) power supply for the production 
of electrolysis, end (2) control and switching equipment for application 
of the electrical power to the electrodes at preselected frequencies and 
pulse durations. 

The negative terminal of the power supply unit was connected to 
the tungsten wires through a relay switch and the positive terminal die 
rectly to the aluminum anode with the dec circuit completed through the 
water, The best results were obtained with a potential of 325 volts 


across its output terminals, 
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The relay switch was activated by @ square weve input from the 
pulse generator at frequencies that were regulated by an ultra-low 
frequency oscillator. 

The operation of the system was controlled manually by a switch 
installed in ths pulsing cireuit between the relay and pulse generator, 
Tho electrolysis circuit was provided with a selector and polerity re- 
versing switch which allowed any combination of the three tungsten wires 
to be switched into the circuit simultansously. A photograph of the 
electrical equipment is given in Figure 3.7. 

For visual observation and photography of the hydregen bubbles 
adequate lighting and a dark background were found to be essential. 

The illumination of the bubbles was provided by two 650 watt 
movie lights, The two light filaments were arranged end to end and 
attached face down to the top of a 12 minch long by 8 inch high by 8 
4nch wide plywood box. The interior of the box was painted with non- 
reflecting bleck paint. A metal partition was installed inside,4 inches 
eway from and parallel to the bottom of the box. To coliimate the light 
1/8 inch wide by 9 inches long slits were cut in the partition and the 
bottem plates of the box. The box was placed on a frame directly above 
the bar supporting the three tungsten wires and adjusted to direct the 
curtain of light from the slits through the plexiglass downward into the 
plane of the hydrogen bubbles, The curtain of light jliuminated approxi-e 
mately one half inch wide strip of bed-material. 

The hydrogen bubbles were photographed through the observation 
window with a 35 millimetre camera equipped with a macro-lens. 

Throughout this experimental investigation a 10 cycle per second 
frequency was employed to pulse the wires, The pulse duration that was 


found to yield best results was approximately 10 milliseconds, To insure 






















| | 
ef? sor] Sogat ovine ofte; CO St SOE ame dod byre ya ier aff arn | 


oy J00¥ eetorompet ga soteveaey eallig 
Liat [Ltese yore = 4 
dodiua a gf fieatin To n0o .Oy ney OM To cabdegéga afk 34a 
cosets} ‘ca yale: wl coded (Rete gufelog a ak cottage aa 
Ot yiualog bit sodoalen « chin Seohewnd aw Sheets ieylonieete ae 
ee crus eds to nolianthne gut Rowlie dotdw dedive gain 
af? Yo te “A teaser Sinoalo ent ojnt bedotive od et 
| Nat pewty a2 Seempbape, Snail 
( ht! nottoreerds Daeaivy egf 
$3 [jaw Dobe goad ald. & Gee geod et eiaepee 
( ee ecloelvd af? te medd Satmunl le edt ee - 


uiewe> ened thedeliey jail =) ef? eddy huaboie 

ol devia Sf 9 Ye qed a? ef oweb ena? Sateen 

non afiw bot | ae xo¢ off So talyeder et? goad booweia obiw doat 

be dhe ok ‘Bul . 0 ew: Solem A deteq Jteald gabteeking! 

ule td ofan? Ile . et af? Ye d of? of. infiwog bee aot ems 
edt bas gier?itueg o4¢ 1? tus otee chtty emak endtoat a A oboe some OME 
eveds yYitemtit exe) 3 ap beoala cay xed ed? ted add to sed at sama 
eat Jowtlt. of bedantine a0 tothe ratagaed eect. od? subteoquee =< 
edd ede8 brawl aea(airele of? dyeoed? Ohta ed wort Stig tf . “— : 

sta ReeadaC? Sight Ye atatele ot spalddvd avgorbyl o 1 aa 


ewmedithrete 9 | 
fen sie + a eo MS ey >a 














Figure 3.8. 





Photograph of Pulsing Prequency Verification 
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that tho applied frequency was constant, pulsation frequency checks were 
made using @n osciillescope, A typical oscillescope record of pulse du- 
ration, its shape and frequency is given in Figure 3.8. This procedure 
confirmed the ability of the electronic equipment to produce constant 
voltage pulses at equal time intervals and of equsl durations. 

The quélity of hydrogen bubbles depended very much on the 
cleanliness of the tungsten wires, Before every experinentel run all 
wires were cleaned using cotton swabs on 6 inch long wooden reds 4nd 
wetted first in Legroine and then in acetone, Depending on the condition 
of the water in the flume, cleaning of the wires was often necessary also 
during test runs, It was found thet the wires could be cleaned with the 
cotton swabs up to one half of an inch from the bed without disturbing 
the bed during a test run. Improvement in the quality of the bubbles 
could also be obtained by pulsing the slectrodes for shert periods of 
time under reversed polarity conditions. 

Often, however, the quality of tho bubbles produced by 4 wire 
varied for no apparent reasen over short periods of time. Such behaviour 
was noted also by Schraub, 2t al. (1965) and has never been explained. 

The life of the tungsten wires was exceptionally long. Only one 
set of 3 wires was used through the entire experimental program. It must 
be noted, however, that no continuous bubble production was ever 
attempted with these wires. 

In the specific application of this techniques to study sediment 
transport phenomen’s additional problems were encountered, ‘Since the 
electrode wires were stretched vertically and passed through the full 
depth of sedimont, hydrogen paspushews ts produced along the length of the 


wire within the sediment and could be seen erupting upward along the wire. 
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To reduce partially the effect of the bubbles generated below the sedi- 
nent surface upon the bubbles in the flow above the lower portions of 
each of the wives were coated with a thin layer of insulation. The 
eruption of the bubbles from below the sediment interface did not occur 
4vmodieately efter the first pulse but was delayed by the time necessary 
for the bubbles to rise into the flow zone. Generally four or five rows 
of bubbles could be pulsed and photographed before eruption effects 
became noted. 

A major difficulty was experienced in obtaining satisfactory 
photegraphic imeges of the rows of hydrogen bubbles near the sediment 
boundary, It was caused by combination of slight bed-wave slope varie 
ations in the lateral direction end reflection of light from the white 
sediment material, On numerous photographs the near-boundary bubbles 
were photographed against a white background of reflected light and 
eould not be identified in the photographs. 

The accuracy of the hydrogen bubble method of velecity measuree 
ment has been investigated by Schraub, et al. (1965). Errors in velocity 
estimation were found to be caused by a number of factors, namely: 

(1) measuroment of displacement between rows of bubbles, (2) averaging 
effects in changing velocity field, (3) displacement of bubbles out of 
the plane of the generating wires, (4) bubble response to fluctuations, 
(5) resolution problems due to finite bubble size and due to finite 
averaging intervals, (6) bubble rise velocity, and (7) velocity defect 
behind bubblo-generating wire. According to the findings of these in- 
vestigators, the total uncertainty interval (estimated confidence inter- 
val, about the expected value of velocity at some specified odds) in 
estimating the true seperation between two adjacent bubbles or bubble 


rows is controlled by the errors of measurenent of the separation dis.. 
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tance, by the averaging uncertainties inherent in evaluating Eulerian 
velocitiss by calculation of Lagrangien time-average velocities, and by 
the displacenent of bubbles out of the plans of the generating wire. 
Other factors were found to introduce only minor uncertainties, In an 
wesarre cited by Schraub, et al. (1965) an uncertainty of 6.2 per cent 
was estimated when the single frame method of data reduction was used 
end bubble displacement was measured 40 diameters downstream from the 
generating wire. 

Although the same single frame method of data reduction was 
employed in the present oxperinents, only the relative varietions of 
velocity were of interest. Hence, nsither the averaging effects nor the 
displacement of the bubbles from the plane of the wire were considered to 
be of sufficient importance and @ 5 psr cent accuracy in the results of 
velocity uedeuremornte was considered obtainable in the present investi. 
gation. 

3.4.4 Pressure Measurement 

fhe pressure measurement system consisted of five components: 

(1) pressure sensor devices, (2) menifold, (3) reference pressure source, 
(4%) two pressure transducers, and (5) recording apparatus. 

The pressure sensors were made of 0,120 inch outside diameter 
(0,012 inch wall thickness) stainless steel tubes. To prevent sediment 
particles from entering and clogging the entire pressure measurement 
system a fine stainless steel screen (80 mesh) was pleced over one end of 
the tube and was soldered to its outside wall, Five or six 0.75 millimetre 
holes were punched in the screen to provide adequate size pressure port 
openings, The other end of the stainless steel tubes was threaded to 


fit vertically into fittings previously installed in the bottom plate 
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of the flume. A total of twelve pressure sensing tubes ranging in height 
from 1/2 to 3-3/4 inches were installed in the flume and were connected 
to a twin manifold system with 1/8 inch inside diameter flexible tubing. 

The purpose of two manifolds was to permit sequential sampling of 
pressures at six sensor positions with a single pressure transducer. In 
addition to the six flume pressure source connections, each manifold was 
equipped also with connecting ports for a pressure transducer, differ- 
‘ential refersnce pressure source, a calibration pressure source or 
system flushing reservoir, and an interconnecting line to the second 
manifold. To monitor pressure from any sensor probe, all redundant 
pressure lines to the manifold were clamped tight with simple laboratory 
clamps, The manifold assembly is illustrated in Figure 3.9, 

A constant pressure source was provided by a free surface water 
reservoir located behind the flume and attached to the laboratory build- 
ing wall. The cylindrical reservoir tank, six inches in diameter and 
eighteen inches high, had five taps each equipped with a valve and a 
fitting for flexible tubing through which the reference pressure could be 
imposed on the two transducers or upon the manifold, As initial obser- 
vations showed that the expected maximum range of pressure variation 
under experimented conditions would not exceed 0.02 inches of water 
(differential pressure) variable reluctance differential pressure trans- 
ducers (PACE Model P90D), rated at plus or minus one inch of water differ. - 
ential, were selected for pressure measurement. A variable reluctance 
transducer has a magnetically permeable stainless steel diaphragm and 
two pressure chambers--one on each side of the diaphragm, The diaphragm 
is clamped between two blocks and deflects when the applied pressures in 
the chambers are unequal, The deflection of the diaphragm changes the 


gap of the magnetic flux path on the system of core and coil assemblies 
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Figure 3.9. 
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that are embedded im each block, The magnetic reluctance varies with 
the gep and produces changes in the inductance of the twe coils. This 
change is converted through spscial circuits of a transducer indicator 
device into veltage representing pressure differences, 

The time variation ef the output voltage from the transducer 
indicaters, which represent differential pressures, were recorded by two 
XwY pen recorders on charts 10 by 15 inches in size, 

Eech of the two separate sensor-.wanifold-transducer-recorder 
systems was thoroughly evaluated to disclese any system anomalies during 
experiments which might have rendered the pressures measurements invelid, 

As illustrated in Figure 3.10, under static calibration conditions 
the transducers exhibited completely linear transfer function with hardly 
@ trace of hysteresis, Gags resolution was found to be limited only by 
the noise in the electrical systems. Using X<Y plotter pen displacement 
of ons inch to represent 0.01 inches of water differential, static differ. 
ential pressure measurements of 0,001 inches of water were consistently 
reproduced. 

An investigation was carried sut also to determine if the abili- 
ty of a pressure sensor to sense pressures below tha water-bed interface 
was affected by the sediment particles surrounding the pressure ports. 
The test was made by recording simultaneously the time variations in 
pressure due to én artificially induced surface wave in the flume and 
sensed by two probes locatsd at the surface of the sediment bed and 
buried three inches below. In all tests, the pressure records from the 
two positions wore found to be identical. 

3.4.5 Observation Area 

Ail instrument probes were installed in an observation area 


approximately 50 feet from the upstream end of the flume. 
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Figure 3.10, 
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To improve visual mesa ntsc acts observations in the flume, 2 
Six feet long plexiglass section of the side of the flume was removed and 
replaced by 3/8 inch thick glass plate, 

All instrument probes were supported by a separate 27 inches long, 
16 inches wide, and one half inch thick plate. It was reinforced by a 
oné inch wide end one half inch thick frame, and was installed at the 
bottom of the flune leaving a space of one half inch between its undersids 
and the flume fleor plate. 

A system of plexiglass bars was installed at the top of the flume 
to support depth and hydrogen bubble wires, and hydrogen bubble illvmine- 
tion lights. . 

The Layout of all probes on the instrument plate is shown in 
Figure 3.11, 

To facilitate photography of the hydrogen bubble all stray light 
and reflections were eliminated from the observation area, Except for a 
feet long glass observation window, both sides of the fluns were covered 
on the outside with non-reflecting black paper and a black curtain was 
installed behind the observer's position to eliminate reflection of the 
laboratory lights from the cbservation window glass. 

The arrangement of equipment, and the working space at the obsor. 


vation window are shown in the photograph given in Figure 3.12, 
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CHAPTER IV 
BED-FORM DEVELOPMENT EXPERIMENTS 
"el Ob 


Sectives and Methods 





Perhaps the most critical part of this investigation was the dee 
velopment of two-dimensional ae formed, equilibrium bed-forms for 
observations and measurements, The difficulties arose because the 
ayteeection processes between the flow, sediments, and the side walls of 
ene flume naturally preduced irregular, three-dimensional bed features. 
Becéuse of this a qualitative study cf bed-form initiation, growth and 
subsequent distortion processes on @ level sediment bed was carried out 
with the following two objectives: (1) to disclose in qualitative terms 
the basic causes for the appearance of irreguler threeedimensional 
bed-features in the flume, and (2) to develop a method by which trains 
of two-dimensional bedforms could be maintained in the flume for 
measurement purposes. 

ALL observations of bed-form formation and distortion processes 
were made in the experimental flume. For each set of observations the 
bed was levelled flat and the flow was established in accordance with 
procedures described in sections 5,2 and 5.4. The experiments were 
carried out over a range of flow conditions, Initial (flat bed) average 
velocities of flow ranged from 2.0 to 3.5 inches per second, and the flow 
width to depth ratios varied from 4.6 to 14.0. 

Observations of boundary geometry charactoristics and of flow 


patterns at different stages of bed-form development were made mostly 
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visually, Injection of dye with a hypodermic needle was employed tc trace 
the structures of flow ebove the bed, 

Attempts to photograph and to video-record the various processes 
were only partially successful, primarily because attempts to photograph 
the initiation process of bed-form development failed. It was found that 
even oxtrems side-lighting angles of the bed failed to produce shadows 


for outlining low bed relief on white coloured sediment material. 


4,2 Natural Formation of Bed-Forms 





A sequence of three discernible physical processes was identified 
in the formstion of bed-forms cn a level bed: (1) initiation, (2) transe 
verse elongation end downstream transletion of crests, and (3) bed-forn 
train formetion. 

After the particle threshold conditions were surpassed, & general 
but spatially discontinuous movement of groups or "patches" of closely 
spaced bed particles was observed. The "patches" varied in size from a 
few particles to several square inches in area, and appeared to be one 
particle diameter thick, Particles were entrained from seemingly random 
locations on the bed and were transported along the bed by rolling or 
sliding in steps of random lengths. On an initially horizontal sediment 
bed, the transport of particles in "patches" took place in the middle 80 
to 85 per cent portion of the flume. 

The natural triggering for the start of a bed-form on a flat bed 
originated from a sudden stoppage on the surface of the bed of one or 
more particles in the lead of a moving "patch", Repeated visual obser- 
vations indicated that immediately prior te stopping, these particles 
broke contact with the bed, and lifted slightly above it. Particle de- 


celeration appeared to commence at the instant their contact with the bed 
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was pescene Their path, as observed from the side and illustrated in 
Figure 4,1 (a) shows particle touchdown at a relatively steep angle to 
the horizontal. It was estimated that the aistanca between the positions 
of uplift and deposition was in the order from 5 to 10 particle diameters. 
No abnormalities on the stationary bed could be detected visually at 
locations of particle uplift or deposition. The formation of the first 
bed-form on the level bed commenced immediately after the eiceerins 
particles came to rest by a build-up of a single layer of particles above 
the mean bed. The downstream particles of the layer aligned themselves 
perpendicularly to the direction of the flow and formed a short length 
crest, 

From the instant of initiation the bed-form started moving down- 
stream and, simultaneously began elongating its crest toward both side 
walls of the flume. Both processes were the result of entrapment of all 
newly arriving particles from upstream at the crest. Particles deposited 
in the lee of the crest advanced the crest downstream. Particles trapped 
near the ends of the crest extended the length of ts crest across the 
flume. 

Particles, arriving at the existing crest, were observed to rise 
slightly pan then come to rest in the lee of the crest. As illustrated 
in Figure 4,1 (b) the paths of particles contributing to crest elongation 
were observed to be similar to the paths taken by the initial triggering 
particles, 

The supply of material for the initial development of a bed-form 
was provided from the same "patch" of particles which contributed the 
triggering particles. However, not all particles of this "patch" con- 
tributed to the bed-form development. Some were deposited upstream of 
the triggering particles, while others, if too distant in the transverse 
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Figure 4,1, Path of Particles 
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direction of the flume, moved past the triggering particles and were de- 
posited on the bed downstream. This suggested that the momentum of 
individual particles in a "patch" was unequal, and that the entrapment 
zone across the flume extended only a limited distance beyond the ends 
of the crest. 

All initially formed crests were straight end clongated toward 
both walls of the flume at velocities which were estimated to be at least 
one order of magnitude greater than the downstream velocities of migration 
of the bed-fern. 

The formation of & bed-form on a level bed triggered the develop. 
ment of a whole system or of a train of new bed-forms, After the trigger. 
ing bed-form had developed a short length of crest, a new crest was ob 
served to form at &@ distance from one to two inches downstream. After a 
short time delay a similar process was repeated again downstream of the 
new bedeform and so om until five or six bedforms appeared on the level 
surface of the sediment bed. 

The process of formation of individual bed-forms in a train 
differed from the process of formation of the initial crest by the source 
of sediment material supplied to each crest. The triggering crest was 
supplied by particles entrained from the level sediment bed from the 
general area upstream of the crest. In a train all bed-forms derived 
their material from within the area between two adjacent crests, 

The processes of formation and growth of a typical bed-form in 
a train of bed-forms started with the development of the immediate up- 
stream bed.form, which for convenience is henceforth referred to as 
bedeform A, These processes aro illustrated schematically in Figure 4.2. 
Initially bedeform A consisted of a single layer of particles above the 


original level bed. The downstream particles of this layer were aligned 
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across the flumo and formed the crest. 

Immediately downstrean of the crest on the original level bed 
a zone of quiescence was established from which particles were not ene 
trained, However, further downstream, approximately 3/4 inches from the 
crest, thers appeared @ zone of enh&énced particle entrainment from which 
particles were scoured below the original mean bed elevation and supplied 
the material for the formation of bed-form B downstream, 

The depth of scour of particles below the original mean bed 
level downstream of bed-form A appeared to be depsndent upon two factors: 
(1) the local intensity of entrainment stresses, and (2) the period of 


time during which an area of the bed was exposed to entrainment stresses 


paration had occurred at the crest of bed-form A even when the crest was 
only one particle diameter high and that the intensity of entrainment 
stresses was increased in the vicinity of reattachnsnt of the separation 
streamline on the back of bed-form B. The peried of tims during which a 
local area on the back of bed-form B was exposed to the enhanced entrain- 
ment stresses was controlled by the velocity of downstream migration of 
bed-form A which in turn depended upon the supply of material to its 
erest from upstream end the height of the crest, 

Observations indicated that during the early stages of bed-form 
development sufficient momentum wes imparted upon all entrained particles 
for them to reach the crest without intermediate deposition and rest, 
Since every particle arriving at a crest was trapped in its lee and 
advanced the crest &@ distances equal to its diameter, the velocity of a 
bed«form was maximum immediately after its formation when it was one 


particle diameter high and migrated on the surface of the original bed. 
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The time of exposure of an area to the enhanced entrainment stresses 
downstream of @ crest was therefore minimum and the scour of bed-material 
from the reattachment area in the initial stages of bed-form development, 
appeared to be limited to a single layer of particles only. 

As the bed-form A arrived into an area from which a layer of 
material had been scoured, its height instantaneously doubled. However, 
its velocity of migration downstream was reduced by half, because the 
supply of material te its crest was continuous and independent of the cone 
ditions at the crest. The increase in the height of bed-form A shifted 
the separation streamline reattachment position further downstream. The 
reduction of its velocity increased the time during which an area in the 
reattachment zone was exposed to the entrainment stresses, And although it 
was not possible to determine conclusively the effect of the sudden change 
in the height of the bedewave upon the intensity of the entrainment 
stresses, the combined intensity-sxrosure effect appeared to enhance the 
initial scouring action, Additional layers of particles were entrained 
from below the criginal bed level and were transported in single steps to 
the crest of bed.-form B. Inmediately, the distance between the two adjacent 
erests started to increase because the velocity of the crest of bed-form B 
remained unchanged and was higher than the reduced velocity of bed-form A, 
However, the relatively high intensity of the initial scouring activity, 
which followed the increase of the height of bed-form A, diminished gradu- 
ally as the process progressed. Newly entrained particles wore not im- 
parted with sufficient momentum to reach the crest of bed-form B in one 
step and were deposited on the back of the bed-form B, Because the re- 
duction in the initial momentum imparted upon the particles and the in- 
crease in the distance between the position of entrainment and the crest 


ef bed-form B were gradual, the deposition of the momentum-deficient 
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particles above the mean bed level commenced at the crest and progressed 
in the upstream direction. The growth in height of the bed-forms thoree 
fore was & result of scour and deposition of particles below and above 
the mean bed level respectively. 

. A train of partially developed bed-forns, photographed 85 minutes 
after the initiation of the first bed-form, is illustrated in Figure 4,3, 
Approximate distances (in inches) of individual bed-form crests in this 
train of bed-forms and their estimated heights (in inches), measured at 
40 and 85 minutes after the initiation, aro given in Table 4,1, The dis» 
tances x were measured from the initial position of the first bed-form 
along the longitudinal centreline of the flume, 

TABLE 4,1 


OBSERVATIONS OF BED=FORM POSITION (x) 
AND HEIGHT (h) VARIATION 
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The above photograph and data illustrate the basic aspects of bed- 
form train formation. During the first 40 minutes after initiation four 
additional bed-forms developed in the flume. Although in time the heights 
and wave lengths of individual bed-forms increased, the rates of increase 


diminished. The velocities of migration of the bed-forms in the downstream 
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direction also decreased, The velocity of a bed-form was, however, 
always greater than the velocity of the adjacent bed-form upstream. 

Although the mechanisms previously described relate to the initial 
stages of bed-form train development they are nevertheless characteristic 
throughout the entire bed-form train development process, The initial 
mechanisms, however, exaggerate some ef the causes and effects and pro@} 
vide better insight into tha entire process, 

Considerations ef the various processes show that the entire bed- 
form system tends to assymptotically approach an equilibrium of constant 
bed-form heights, wave-lengths, and velocity of migration. When a newly 
formed crest one particle diameter high doubles its height and reduces 
its velocity by half, entrainment stresses, their location and area 
exposure tins are suddently and drastically altered in relation to 
previous conditions cownstrean of the crest. However, when a hed-forn 
is several particle layers high, and the change in the height is by one 
particle diameter only, the relative effects upon the particle entraine 
ment parameters and, hence, upon the geometry of the downstream bed.form 


are much less sévers, 


nts 


4.3 BedeForn Distortion Expsrins 





4.3.1 Introduction and Classification 

The formation of irregular three-dimensional bed-features on an 
initially level sediment bed was observed to progress through three 
stages of development, Each stage was distinguished by a distinctive 
bed-geometry, namely: (1) two-dimensional bed-forms, (2) symmetrical 
three-dimensional bed=waves, and (3) irregular three-dimensional bed-= 
features, 


In the bed-form distortion experiments, the development of bed- 
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forms was induced artificially by the wire method described in Section 
5.4. No discernible differences were noted between the development of 

| bed-forms downstream from an artificially induced ripple and naturally 
triggered bed-forms. Up to seven two-dimensional bed-forms were ob- 
served to form downstream of the initiating ripple before the appearance 
of three-dimensional effects, 

A simple investigation of flow above a two-dimensional train of 
bed-forms was carried out in conjunction with the distortion experiments, 
Coloured dye was injected into the flow at various locations above the 
bed. It was found that in a 1/4 to 1/2 inch wide zone adjacent to the 
side walls of the flume the flow was laminar, In the remainder of the 
flow, completely turbulent conditions were indicated by a rapid 
diffusion of dye into the flow. 

4.3.2 etrical Three-Dimensional Bed-Features 

The intermediate stage of bed-form development was dis- 
tinguished by three-dimensional bed features which were Symmetrical 
about the centreline of the flume, Although generally, distortion 
of bed-forms may be caused by numerous other factors, in these experi- 
ments the distortion was caused by side-wall effects, 

On a fiat bed, the average flux of particles was uniform in the 
central, 80 to 85 per cent of the flume only, and reduced gradually to 
zero at the side walls, Because of the reduced particle supply near 
the walls, the central portion of a bed-form crest translated downstream 

at a velocity greater than its ends, the crest assumed a flattened 
U-shape and merged into the flat bed near the walls as its height 
gradually decreased to zero. 


Downstream of a bed-form, entrainment of particles was always 
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enhanced slong the separation streamline reattachment position. From 
the reattachnont area imnediately downstream of the curved portions of 
the crest, sediment particles were transported toward the walls and the 
centre of the flume and were deposited forming longitudinal ridges with 
depressions between them, Together with portions of the two adjacent 
crests, these ridges formed pairs of characteristic "bowl" features. 

A typical intermediate stage bed geometry is illustrated in 
Figure 4,3. The flat bed width to depth ratio was 4.6, and the average 
velocity, computed from discharge meter reading, was 2.9 inches per 
second. The photograph, taken 85 minutes after the upstream bed-forn 
was artificially initiated, illustrates symmetrical, three-dimensional 
bed«features, The crest of the first bed-form is straight and perpendi. 
cular to the flow in the central portion of the flume only. It curves 
backwards from approximately 4-1/2 inches from each wall, Two "bowl" 
features are seen on the back of the second bed-form and have affected 
the crest of the downstream ripple. The material, excavated from within 
the "bowls," and transported to the crest along the inner ridges, has 
lecally advanced the crest and has produced protrusions on each side of 
the centreline downstream of the crest. The third and all subsequent 
bedforms sean in the phetograph, have also developed thrse-dimensional 
features, Their crests have been retarded near the walls and near the 
centreline of the flume and the pairs of inner longitudinal ridges, such 
as are seen on the back of the second bed-form, have joined together and 
have formed the high ground near the centre of the flume. 

The photograph shown in Figure 4.4 illustrates bed-features when 
the initial width to mean depth ratio in the flume was 9.0 and the average 
flow velocity was approximately 3.2 inches per second. The "bowls" seen 


on the back of the second bed-form have been elongated, thoir inner 
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Photograph of Intermediate Stage Bed-Form 
Development 
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ridges have approached each other, and the crest line, while still 
Symmetrical, has acquired reaches of substantial curvaturs, 

The structures of flow, associated with the intermediate stage 
of the bed-form distortion processes, were investigated by tracing the 
movement of dye injected at various locations above the bed, and by 
observing visually the movement of sediment particles along the bed, 
Some observations wers recorded en video-taps for repeated viewing, The 
results of those observations are summarized in Figure 4,5, 

Excluding the laminar flow conditions which existed near the 
Side walls of the flume, it was found that two distinct zones character. 
iced the flow in the intermediate distortion stage: (1) an upper zono, 
and (2) a lower zone, The flow in the upper Zone was continuous over 
all bedefeatures, and was characterized by @n absence of macroescale 
flow structures. No permanent transverse velocity components in this 
zone couid bo identified. The flow in the lower zone was discontinuous 
and extended downstream from the crest to 4 position on the back of the 
downstream bed-form. I+ was bounded by the uppsr zone flow above and 
the sediment bed below. However, its primary distinguishing characteris- 
tics were the macro-scale flow structures--the rollevs and the spiral 
flew cells. 

The rollers were present in the lee of straight reaches of 
bed-form crests, Their presence was indicated by &@ small net displace. 
ment of sediment particles in the upstream direction toward the toes of 
bedeforms. Instantaneous diffusion of injected dye indicatsd that the 
intensity of turbulence in the roller was relatively very high. In the 
&rea downstream of the curved segment of a crest, close to the sediment 
bed and near the toe of the bed-forn, a permanent drift of fluid toward 


the centre of the flume was indicated by dye injection. The velocity of 
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this flow was very slew and flow conditions were considerably below 
particle threshold value. The dye diffused into the main stream only 
efter it moved over the line extending upstream from the inner longitu- 
dinal. ridge, 

The axes of the spiral flow cells were eéligned in the general 
Cownstream direction. Whenever a pair of "bowl" features had Geveloped 
on the back of a ripple, a total of six Spiral cells, three of each "bowl", 
could be identified. Two of them occurred within the "bowl", and the 
third was found to exist outside and édjacent to the inner longitudinal 
ridge. The spiral cell located within a "bowl" and adjacent to the inner 
Longitudinal ridge was found to be more intense than the other two and 
could be more easily demonstrated by dye injection, It was designated as 
the primary spiral. For an observer looking downstream, the directions 
of rotation of the primary spirals were counterclockwise for the left 
"bowl" and clockwise for the right "bowl", The direction of rotation of 
the two secondary cells, associated with each "bowl", were of opposite 
senses to those of its corresponding primary spirals, The vortex-like 
features of the flow disparsed gradually into the upper flow zone before 
reaching the crest of the downstream bedeforn, 

The general directions ef fluid motion in the macro-scale 
structures strongly suggested that a maton axis of rotation joins 
ths rollers with the pairs of the spiral cells in the individual "bowls" 
or in the central portion of the flume downstream of the straight reaches 
ef the crests, Attempts to demonstrate the continuity of the rotation 
éxis, however, were not successful because of rapid diffusion of injected 
dye due to severe intensity of turbulence in the zones where the curvature 
of the axis of fluid rotation was most pronounced, 


The longitudinal ridges on the bed were built up by particles 
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delivered by the transverse flow components near the bed. As time pro- 
gressed, the height of theses ridges, particularly the inner ores, in- 
creased and eventually appsared to reach a maximum. From thet time, all 
particles delivered by the spiral flow components to the crests of the 
Praga were immediately swapt downstream. The flux of particles was 
therefore more intense along the ridges, than in adjacent locations and 
affected the aerial geometry of the downstream bed-form by non-uniform deo 
position downstream of the crest, 

4,303 Bed-Form Distortion 

The onset of loss cf symuetry in the boundary forms signified the 


appearance of the final stage processes responsible for distortion of the 


The distortion processes usually originated in the vicinity of the 
third or fourth ripple in 4& train, locally, their effects on the bed 
forms were almost instantansous, and they progressed downstream effecting 
established bed-forms as well as the flat bed, The rate at which these 
effects propsgated downstream and established irresular boundary forms on 
the bed was estimated to be at least one order of magnitude faster than 
the development of ripples in the initial and intermediate stages, 

The triggering mechanisms, which initiated the distortion processes 
originated generally at locations where a slight asymmetry of boundary 
features had developed, These were most severe in the vicinity of close- 
ly spaced or joined inner ridges, or in the area downstream of large 
curvature of crests such as existed downstream of crest protrusions. 

An investigation of flew structures above the distorted boundary 
was inconclusive, However, transverse flow velocity components were de- 
tected throuzhout the entire depth of flow, indicating thatthe two distinct 


depth zones of flow, associated with the intermediete stage, did not exist. 
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Visual observations of the movement of particles, Suggested also thst 
transport was not uniform across the lames Some particles were entrained 
and moved on the backs of the bed waves, Other particles moved in higher 
concentrations along selected routes between the irregular three- 
dimensional bed-forns,. | 

In addition to the distortion mechanisms described in previous 
paragraphs, symmetrical bedforms were cften disturbed by other causes, 
Typical of these were the Veshaped ripples appearing on the backs of 
existing bed.forms, They originated at slight depressions on the bed and 
developed two straight crests which in plan represented the spreading 
arms of aV, At high velocities of flow this type of disturbance often 
resulted in a division of the criginal bed-form into two or more separate 
bed=festures,. 

4.3.4 Analysis of Observations 

The main reasons for the distortion of bed-ferns into irregular 
thres-dimensional patterns were traced to the side wall effects upon the 
flow. When, andes generally turbulent flow conditions, the particle 
threshold was reached throughout most of the bed, the aA nS eee layer 
et the plexiglass walls wes laminar and a transition from laminar to- 
turbulent flew was established in a narrow zone adjacent to the wall, In 
this zone of transition at the bed level the rates of entrainment and of 
transport of particles per unit area or per unit width psrpendicular to 
the mean flow reduced from a maximum at the fully turbulent edge to zero 
close to the wall. Ths variable rate of particle transport contributed to 
the noneuniform rate of crest movement Howiatrecn across the flume, and 
affected the vertical growth rates of the bed-form, As a result, the 
crest of a newly developed bed-form was straight and the bed-form was 


of a uniform height in the central portion of the flume only. Near both 
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walls, the crest line curved backwards and the height of the bed-form 
was reduced gradually to sero at or near the wall, 

Because of flow separation all along the crest line, the separation 
streamlines defined a complex, three-dimensional surface in space. The 
mean intersection line of this surface with the sediment bed surface on 
the upstream slope of the adjacent bedeform in the central portion of the 
flume was perpendicular to the mean flow, Near the walls, however, down- 
stream of the curved segments cf the crest, the reattachment line was 
curved and induced flows with transverse components across the flume. 

Although it could not be determined if secondary circulation had 
not been established in the flume prior to the formation of the first 
bed-form, the appearance of spiral flow cells was dirsctly linked with 
the establishment of crestline curvature and the associated transvere flow 
components. 

A local net movement of sediment in both directions from the 
general Bree aan’ area at the reattachment lines was produced by the 
near.boundary flow. From the curved segments of the reattachment line, 
one portion of the entrained particles was transported toward the nearest 
wall, and the other toward the middle of the flume. The ability of this 
flow to sustain particle transport, however, diminished rapidly from the 
entrainment area, 

The initial deposition of particles near the walls appeared to be 
caused by the upward deflection of the neareboundary flow as it approached 
the wall. Particles came to rest because the flow could not sustain 
their movement vertically upward against the gravity. The reasons for 
the deposition of particles transported toward the centre of the flume 
are less clear. At least three causes, Slone or in combination with each 


other, may have effected the transporting capacity of this transverse 
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flow, To maintain the continuity requirement, the near-boundary stream 
could have been deflected upward by the influence of its counterpart 
approaching from the opposite side of the flume. It could also have been 
retarded by the development of an interface with the flow downstream of 
the straight portion of the crest or by boundary friction. 

The development ef secondary spiral flow cells and the longitudi- 
nal bed-features appsared to be interrelated in a reciprocal fashion. 
Although the presence of the spiral cells wes demonstrated only above 
bedefeatures of considerable size, the continuous developnent of the 
Longitudinal ridges from the Boca when curvature of the crest was very 
large ena sted that the secondary cells were also established at that 
time. 

As time progressed, the bedforms increased in size and the 
secondary spiral cells increased in diameter, More particles were en- 
trained by the enhanced erosive action along the reattachment lines and 
were moved toward the ridges by the transverse flow components inherent 
in the spiral cells, In addition, previously deposited particles were 
eroded from the sides of the inner longitudinal ridges and the ridges were 
Slowly displaced toward the middle of the flume. 

The rotating flow of the primary spiral cell separated from the 
boundary as it approached the crest of the inner ridge and a counter cir- 
culation of lesser intensity (the secondary spiral cell) was developed on 
the opposite side. The direction of rotation of the fluid suggested a 
U.shaped vortex tube made up of the secondary spiral at the two inner 
ridges and the roller in the lee of the straight portion of the crest. An 
attempt to establish continuity of the axis of this vortex was, however, 
unsuccessful because of the low intensity of fluid rotation in the roller 
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and the high intensity of turbulence at the expected junction locations 
of the roller and the Spirals. 

During the transitional period of bed-wave development, charac- 
terized by the symmetrical pattern of bed-features, the longitudinal 
extent of the primary and secondary spiral cells appeared to be controlled 
by the upper flow. Judging from the direction of particle motion in the 
primary spiral cell area, the near-boundary flow components were nearly 
perpendicular to the longitudinal axis of the flume near the point of the 
origin of the spiral and nearly parallel to the axis of the flume where 
the spiral motion merged with the upper zone flow. The rotating fluid 
therefore appeared to have been accelerated by shear at its upper bounda- 
ry and was diffused into the upper zone flow before reaching the crest of 
the bed-form. The rotational motion features of the secondary spiral 
cells were similarly affected. 

Although the direction of motion of particles approaching the 
crest of a bed-form was generally downstream across the full width of the 
flume, zones of particle transport concentration were nevertheless estab. 
lished on the backs of bed-forms with symmetrical bed-features. The zone 
of maximum concentration occurred along the inner ridges, Particles en- 
trained in areas where transverse boundary currents were present were 
first delivered to the longitudinal ridges and only then were transported 
downstream, 

The variation of concentration of moving bed particles across the 
flume affected the crest line of the bed-form at which all particles were 
deposited. The high concentrations of particles transported along the 
inner ridges, locally increased the rate of the crest movement downstream. 


The crest line acquired segments of additional curvature and the separation 
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streamlines produced & separation surface of even higher complexity. 

From the observations it was concluded that during the transi- 
tional period a semi-equilibrium state existed between the mutually com. 
patible system of symmetrical bed~features and the secondary flow struc 
tures above, Both, the bed-features and the secondary flow cells, ro- 
tained thelr individual identities in the total system from the outset 
of their formation, but continued to increase in physical size and became 
more sensitive to small disturbances as the processes ae, The 
spatial extent and influence of the secondary flow cells was restricted 
to individual bedforms and appeared to be governed by the flow above. 

The disruption of the symmetrical bed-features and the subsequent 
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small disturbances of the symmetrical geometry of the bed or of the sec- 
ondary flow structures, 

Typical disturbances were generated downstream of a bedeform crest 
which was subject to locally accelerated rates of particle deposition. Be» 
cause of ths modified crest line curvature, the separation surface inter. 
section with the bed produced local transverse currents that presumably 
were not compatible with the existing system of secondary flow cells. A 
different type of disturbance was generated whon a pair of inner ridges 
on the back of a bedewave approached each other and the associated secondary 
spiral cells merged together. In all cases the restructuring of the flow | 
produced immediate effects on the bed downstream, disturbing the symmetri- 
cal bed-features into irregular shapes, 

The observations of the sequence of the events, from the initia. 
tion of bed-forms to their subsequent distortion, suggested that the 


different stages of their development characterized by the symmetrical 
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énd irregular features represented two distinct visible phases (Rlench, 
1969) of bed-forms that generally corresponded to the ripples and dunes 
as described by Simons and Richardson (1966), The imposed flow, in all 
cases, led to the development of dunes, The aatemeiieee stage of sym- 
metrical bed-forms represented a transitional state tending toward an 
equilibrium in the duned range, The distinguishing characteristics bee | 
tween the two phases, based entirely on simple visual observations may 
be summarized as follows: 

(1) Ripples are associated with small bed-wave heights and small 
scale secondary flow cells that =e bounded by a continuous upper flow - 
zone that extends in both directions, perpendicular and parallel to the 
flow across the flume. The effects of tho secondary flow cells are not 
transferred across individual bed-wave boundaries, but through shear are 
dispersed into the upper flow befere reaching the crest of the bounding 
bedewave, 

(2) Dunes are associated with larger bed-form heights, and at 
ieast some of the secondary flow cells extend not only the full depth of 
flow, but extend over a number of adjacent dunes in the direction of flow, 
The continuous secondary cells tend to oscillate fron Side to side in the 
flume, isolating local zones of flow separation from crest of three- 
dimensional bed-waves, And, sediment transport takes place not only on 


the backs of the dunes, but along continuous routes of tho secondary cells, 


4,4 Bed-Form Control Technique 
From the findings described in the previous section it was con» 


eluded that the distortion of sediment bed-features into three-dimensional 
shapes was caused by the local retardation of the flow near the side walls 
-of the flume, Local retardation of flow reduced transport of particles 
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and produced bed-foima crest curvature and cells of secondary flow whose 


stability eventually determined the stability of the total system. 

The development of irregular three-dimensional bed-forms was 
successfully eliminated by preventing the development of excessive crest 
curvatures near the walls of the flume. A manual technique was employed 
to maintain a straight crest across the flume. Whenever the central por- 
tion of a bed-form had advanced a distance of approximately one inch 
ahead of the crest at the walls, a 1/16 inch diameter wire was inserted 
vertically into the sediment, approximately one inch deep, one inch up- 


stream of the crest, and 1/4 inch from the wall. It was then moved 


t > Ae) 3 uv 2 
Slowly downstream dragging adjacent sediment particles forward until a 
4 < q me q ae Cade eae ae - i = ~ 
straight crest line was restored, For best results the spesad of “he wire 


had to be relatively uniform and slow. The resulting bed-forms are 





Figure 4.6, Photegraph of Bed-Forms 
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The effects of wire insertion into the sediment bed upon the 
flew and pressure fields at the instrument probes were investigated 
separately, Pressures were monitored and transducer output was ree 
corded continuously while bed-form crests were straightened by the wire 
method. The results showed that, provided the wires were inserted 
vertically and were moved gently, pressure fields at sensor positions 


remained unaffected. 
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CHAPTER V 
EXPERIMENTAL PROCEDURES 


5.1 Outline 

A description of data acquisition procedures is given in this 
chapter, The presentation is divided into five topics which follow the 
chronological sequence of procedural steps in a typical test run: 

(1) Flume and sediment bed preparation. 

(2) Preliminary instrumentation check and flow depth adjustment. 

(3) Bedewave development method. 

(4) Observation and data acquisition procedures, 


(5) Data reduction and handling. 


5.2 Flume and Sediment Bed Preparation 

52.1 Levelling of Sediment Bed | 

To insure identical initial conditions in the flume in all tests, 
standard flume and sediment bed preparation procedures vere evolved and 
used throughout this investigation. 

One of the preliminary tasks was the preperation of a four inches 
deep horizontal bed of sediment from the upstream end of the flume near 
the headbox to approximately twenty feet downstream from ths observation 
arsa, 

The method for levelling the sediment bed was discovered by chance, 
It was observed that when the water-sedimont mixture in the flume was 


oscillated across the flume at its approximate natural frequency, and then 
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was allowed to dampen naturally, that the surface of sediment particles 
tended to assume a horizontal level, 

To level the sediment bed, the surface of the water in the flume 
was initially adjusted to approximately 1/4 inches above the required 
horizontal bed surface level. This submerged all granular particles and 
provided reference plane from which small variations on the bed surface 
could easily be observed through the sides of the flume, 

To level the sediment bed, a 1/32 inch thick, 18 inches long, and 
10 inches wide aluminum plate was used, The 18 inch axis of the plate 
was aligned horizontally with the longitudinal centreline of the flume, 
and it was inserted vertically into the water-sediment mixture. The plate 
was then oscillated by hand with an amplitude of 2 to 3 inches and was 
gradually advanced along the length of the flume, 

To obtain level bed surface this procedure had to be repeated 
several times. Seldom, if ever, the initial material distribution along 
the flume was uniform, and after initial levelling it was often necessary 
to remove or replenish sediment particles from high or low spots. 

The effectiveness of this method was investigated by further 
lowering of the water level in the flume to expose approximately one half 
of the surface particles, Judging from the accumulations of the granular 
particles above the water line it was concluded that the areal distribution 
of high lying particles was not uniform. There was a tendency for the 
exposed particles to accumulate into areas of irregular shape and size. 
However, with rare exceptions, the variation in depth between the highest 
exposed and low submerged surface particles did not exceed two particle 
diameters, This was determined by filling a few random low spots with wet 


sediment particles in lifts one particle thick until a common level with 
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the highest particles was reached, 

502.2 Flume Filling Procedure 

At the conclusion of the levelling process the depth of water 
above the topmost surface particles was only a small fraction of en inch. 
To avoid excessive flow velocities capable of disturbing the surface 
particles during the flume refilling operation, the first one inch of 
depth above the sediment bed was raised at a very slow rats, usually in 
about 6 hours, The filling rate could then be increased and the flume was 
filled with water to a depth of one inch above the anticipated flow depth 
for the test. 

After filling, tho weter in the flume was left for periods of up 
to 48 hours to deaerate and to warm to laboratory temperature, The mu- 
nicipal water used in the flume had tap temperatures less than 40 dogress 
Fehrenheit which was considerably below the norzal laboratery tempsrature. 
Deasration was essential because air bubbles in the water produced scatter 
of light and prevented observation of hydrogen bubbles used in velocity 
measurements, In addition, there was a tendency for the minute air bubbles 
in water to cling to Sndividual sediment particles and thus modify their 


specific gravity. 


5.3 Pre-test Instrument Check and Flow Depth Adjustment 





5301 Verification of the Measurement Systems 

The pre-experinent pressure measurement system verification had a 
dual purpose, First, it disclosed anomalies in the system, such as the 
presence of entrapped air bubbles or blocked pressure probe ports, Secondly, 
it provided a record of a reference pressure for all pressure measurements, 
essential in evaluating the comparing the experimental pressure data, 


The pressure measurement system tests were carried out using the 
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still surface of water in the flume and an auxiliary free surface reser. 
voir 8s differential pressure reference, Equipment verification tests 
were attempted only after all free surface oscillations and surface waves 
in the flume were completely damped out, Depending upon the severity of 
the initial disturbances, the time period required for the free surface 
oscillations to subside ranged up to six hours, during which time no 
activities were permitted on the experimental floor, 

The usual steps in preparing the pressure measurement equipment 
for an experiment were as follows: 

(1) The electronic equipment was switched on and permitted to 
warm up for 1-1/2 hours. Praliminary investigations had disclosed that 
this period of time was necessary to stabilize a very slow, exponentially 
decressing electrical drift in one of the two pressure transducer systems. 

(2) Zero differential pressures were imposed en the prassure trans. 
ducers by interconnecting both pressure chanbers of each transducer to the 
same reference pressure source and transducer output voitages were then 
adjusted to zero, if necessary, 

(3) The positive pressure chambers of the transducers were cone 
nected, through the manifolds, to any two of the twelve pressure probes in 
the flume, and by using the established zero differential pressure as refer. 
ence, water level in the auxiliary free surface reservoir was equalized with 
the water level in the flume. 

(4) Pressure transducer output recorders were started and each of 
the pressure probes in the flume was then sequentially connected, through 
a manifold, to a transducer and a one minute sample of its output under 
the &pplied zero differential pressure was recorded, In addition an out- 


put record was obtained with both chambers of the transducers connected to 
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the auxiliary free surface wank. This procedure was repeated three times 
to ascertain the consistency in individual probe outputs. 

A differential pressure of 0,002 inches of water above or below 
the mean output of the six pressure probes used with each of the trans- 
ducers was arbitrarily selected as the standard indicating no anomalies 
in the individusl probe systems, provided that consistency in this output 
from repeated messurements was maintained, If the output from any pressure 
probe was found to bs outside of this spscified limit, or if variation was 
detected between repeated measurements, the pressure lines to the probe 
were flushed with slightly soapy water threugh the special ports in the 
menifolds provided for in-place transducer calibration, Usually, the 
eauses of these small anomalies were tiny air bubbles attached to the 
interior of the pressure lines or at the junctions of the flexible tubing 
and metal fittings. In most cases the remedial flushing procedure effece 
tively climinated the problem. 

During the initial stages of the experimental program, gage cali- 
brations wore checked before each test was run, This pre-test step was 
omitted as the work progressed in favour of calibration verification 
between test runs, The pressure measurement systems were observed to 
maintain constant response cheracteristics, and the checks served pri- 
marily to insure that the proper range settings in the electronic circuits 
were not accidentally changed. 

After ths conclusion of pressure system verification, flow was 
established in the flume by starting the pump motor and slowly opening the 
flow contrel valve until a discharge required for the test was attained. 

Hydrogen bubble apparatus was switched on, and after a short 
warm-up period, the wires wers pulsed and the production of hydrogen 


bubbies at each of the three wires separately and in the four different 
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combinations with each other, was observed visually to insure proper 
equipment operation. 

503-2 Water Level Control 

_ Preliminary observations indicated that because of evaporation 
from the free surface, leskaze from the flume and withdrawal of cooling 
water for pump bearings small amounts of water were lost from the recircue 
lating system, Although these water losses did not exceed 0.25 eubic ine 
ches per second the corresponding continuous lowering of the mean water 
surface from the flume affected all pressure measurenonts since they were 
referred to a constant pressure source, 

Because the rates of leakage under typical test conditions were 
roan to be constant for periods up to six hours, to maintain a constant - 
mean surface level in the flume an equivalent amount of water was replaced 
inte the system, 

To determine the rate at which water had to be replaced, the flume 
was filled one inch deeper than the required experimental depth and a dis. 
charge of water equel to the required experimental discharge was established 
by adjusting the flow control valve, As the flow in the flume under these 
conditions was insufficient to entrain or transport sediment particles, 
and the bed remained level, the piezometric head at any location remained 
invariant except for the drop in the water surface elevation due to leak. 
age, Before a test was run, therefore, pressure was monitored for 30 mi- 
nites at two probe positions in the flume, The corresponding leakage from 
the flume was calculated from the rate of the recorded pressure drop and 
replaced back into the flume from a constant head tank loceted near the 


upstream end of the flume as indicated in Figure 3.1. 
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5.4 Bed-Forn Development 

Initial disturbances of bed particles by surges induced by come 
moncenent of flow in the flume were eliminated by establishing a pre- 
selected expsrimental discharge at a depth which produced flow velocities 
and conditions above the bed considerably below particle threshold values. 
After all possible surging activity in the flume had subsided the water 
level in the flume was slowly lowered by keeping the system drainage valve 
open until the flow velocity surpassed particle threshold velocity and 
movement of bed particles could be observed, 

To start the developnent of bed-forms at any location in the flume 
& 1/16 inch diameter wire, inserted vertically approximately 1/4 inches 
into the sediment bed was drawn from one side of the flume to the other. 
If the depth of the ponetration of the wire into the sediment, and the 
speed with which it was mowed were uniform and steady, the particles dis.- 
placed by the wire and flow in the flume formed a trough-ridge combination 
of a relatively uniform cross-section across the flums, If this proce- 
dure was carried out after particle threshold conditions on the original 
bed had been surpassed, the wire-formed ridge became the crest of a bed- 
form and bed-forms then formed downstream of each other in a@ manner that 
was identical to the naturally initiated bedforms, 

To eliminate the effects of possible changes in the mean depth of 
flow in the observation area because of boundary roughness changes during 
the measurement interval, bed-form development was triggered at four 
equally spaced locations in the eépproach chénnel, This system of bed- 
forms developed quickly and eventually distorted into an irregular pat- 
tern as no attempt was made to control the shape of the bed-forms manually. 


Once fully developed the irregular bed-forms produced an equilibrium bound- 
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ary roughness that on the average remained constant over the period of 
observation, and produced a constant energy gradient. 

A more severe difficulty was encountered in the immediate ap. 
prosch channel upstream of the observation deck because bed-forms arriv~ 
ing from upstream were three-dimensional, To facilitate restoration of 
two-dimensional bed-forms two 1/4 inch opening sereens were placed 
across the flow, approximately six inches from each other, and ten feet 
upstream of the observation deck, The size of wire of these screens was 
approximately 1/16 inches, and the restriction was sufficient to suspend 
Lomporarily all of the bed particles arriving at the screens, The sus- 
pension of bed-material served two functions, It facilitated a more 
and reduced secondary 


wey 


uniform distribution of particles across the flum 
ae effects in the approach channel. Two-dimensional bed-forms 
were restored using a 1/16 inch diameter wire. -However, restored bed- 
forms. were never used for measurement purposes, but served only to main- 
tain a uniform boundary roughness in the channel, 

Subsequent evaluation of pressure records indicated that bed- 
forms developed in the manner deseribed in previous paragraphs produced 
an average constant depth of flow at the observation area, Provided that 
the quantity of water in the system remained constant, local pressure 
variations were therefore attributable to the passage of the bed-wave 
itself. 

To summarize, the procedure of bed-form development for experi- 
mental observations consisted of the following steps: 

(1) After the required flow was established, the water level in 
the flume was gradually lowered by draining, until the chosen flow depth 


was reached, 
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(2) Bed-forms were initiated at a location approximately 8 feet 
upstream of the observation deck and also at four locations in the 
approach channel upstream of the screens, 

(3) The development of bed-forms in the channel area immediately 
upstream of the observation deck was controlled manually throughout their 
growth period, and then until the experiment was completed, 

(4) Irregular bed-weves arriving from upstream were restored at 
wavelength intervels corresponding to the naturally developed wavelengths 
of bed-waves approaching the observation area, 

In all experiments the maintenance of the bed.waves in the approach 


channel was carried out by an assistant cbserver. The maintenance of bed= 


ried out by the principal observer, 


5-5 Observation and Data Acauisition Procedures 

5. 5el Outline 

Discussion of experimental procedure are divided into two general 
parts. In the first, summary outlines of measurement procedures of indi- 
vidual parameters are given separately, Then the sequence of operations 
and procedures integrating all sepsrate parameter observations is given. 

A number of trial experiments were run to determine most efficient 
procedures to obtain adequate quantity of good quality experimental data, 
Experimental procedures were evolved from experience gained and obser- 
vations made in these trial tests. 

5.5.2 Dopth-Wire Observations 

Bed=wave geometry data was obtained by visual cbservation of the 
instantaneous position cf the surface of sediment at each depth-wire. 


The observed readings were voice recorded on magnetic tape at the time 
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of observation. 

‘ The method cf voice recording of visual observation deta was 
found to be most effective. On the average, it took 25 seconds to scan 
the positions of the sediment bed at all nine depth-wires, As in most 
tests, this scan period corresponded to approximately 1/50 to 1/100 of 
the bed.wave period, local depth corrections because of bed-wave eras 
tion downstream and relatively small bed surface slopes were minimal and 
were not required, 

A typical sequence of steps for bed-wavo profile observation was 
as follows: 

(1) Depthewires were illuminated by ehinate on the hydrogen 
bubble light. 3 

(2) The reference position at the upstream depth-wire was de- 
termined from the known position of the free surface of flow. This step 
was necessary because depth-wire scale had no numeric code to indicate the 
integer number of inches from the reference level to the surface of the 
sedinent, 

(3) Timo from the start of experiment was noted and recorded to 
denote the beginning of the period of depth-wire observation. 

(4) Starting with the upstream wire, sediment level positions 
were observed at each depthewire and were voice. recorded on magnetic tape. 
Positions and depths of bed-wave crests or toes occurring between wires 
were estimated from positions of bed levels at adjacent wires. 

(5) Time on the experimental clock was observed and recorded at 
the conclusion of the scan of depth-wires,. 

5.5.3 Velocity Observations 

In a typical test, all electronic equipment for the production of 


hydrogen bubbles was switched on before the test run was started except 
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for the power unit, which, for safety reasons, was switched on only a 
short period before the wires were pulsed, 

The typical procedure for velocity measurement with the hydrogen 
bubble technique is outlined as follows: 

(1) The power unit and illumination lights were switched on and 
the selected wires were pulsed to observe visually through the side glass 
the quality of bubbles, This procedural step indicated if cleaning of the 
hydrogen bubbles was necessary, This often eliminated a wire cleaning 
process reducing the possibility of local bed-wave disturbance, It also 
indicated if any of the wires produced poor quality bubbles that were 
difficult to observe and would not be recorded on the photographic film. 

(2) If wire cleaning was not required, this step was omitted. 
Otherwise, the dirty wires were first sub jected, for 15 to 20 seconds, to 
reversed polarity pulses, followed by a normal pulsation period of approxi- 
mately the same duration. Experience showed that the reversal of polarity 
was most often sufficient to clean the wires and to enhance the quality 
of the hydrogen bubbles, The reasons for this behaviour, already dis- 
cused by Schraub, et al. (1965), have not been explained and remain une 
known. If further cleaning was required, it was carried out by the method 
described previously in Chapter III with the power unit switches off. 

(3) The camera position was adjusted to cover the field of the 
hydrogen bubbles, The lens was focussed on the hydrogen bubble wires. 

The camera height was adjusted to be at approximately the same elevation 
as the sediment level, This reduced or eliminated the white background 
area, due to illumination of the sediment surface by the light from above, 
against which hydrogen bubbles could not be seen. 


(4) The wire or wires were then pulsed by manually operating the 
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pulsing switch and & photograph of the rows of bubbles was taken by ma- 
nually depressing thecamera shutter before the first row of bubbles pro- 
duced at a wire had reached the position of the immediate hydrogen bubble 
wire downstream, This avoided overlapping of the rows of bubbles origi- 
nating at the two wires, 

(5) Tima from the beginning of experiment, film frame number, and 
the wires pulsed were recorded, 

5.5.4 Pressure Measurenont. 

As described in Chapter III, the instrumentation for pressure 
measurement consisted of two differential pressure transducers, each of 
which could be connected through a manifold to any one of six prossure 
sensors fy the flume. Reference neee. was provided by a constant wa- 
ter level tank which was connected through separate lines directly to the 
negative sides of the two pressure transducers, Transducer output vol- 
tages were recoried in anslog form by two X.Y State See, the X.axis roe 
prosenting time, and Yoaxis the transducer output voltage, 

Because of very small absolute magnitudes of pressure changes in 
the light weight sediment bed, adequate resolution ef transducer output 
could be obtained only at high sensitivity settings of the X-Y recorders, 
This restricted the range of input voltages into the recorders + 0.02 
volts DC, and in effect, restricted the operational differential pres- 
sures between the flume and the referance pressure to the near zero range, 

Pressure measurements in the sediment bed were obtained in two 
ways. In the first, which was used in the majority of experiments, the 
pressures at each sensor position in the flume were sampled for short pe- 
riods of time in @ predetermined or a random sequence. In the second way, 
either one or both of the transducers were connected to a single sensor 


probe in the flume for the duration of the experimental run or for a 
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complete bed-wave passage time period over the sensor, 

The pressure data acquisition procedures included also steps to 
protect the recorders from possible mechanical damage due to excessive 
input voltages. It was found that removal or replacement of clamps on 
the pressure lines produced sudden line pressures and consequently trans- 
ducer output voltages considerably in excess of the available recorder 
pen travel range under the operating sensitivity settings, To protect 
the recorders, a switch was installed in their input circuit, The vol- 
tage output into the recorders was switched off manually whenever differ- 
ent sensor pressure lines were opened to the manifold. 

The pressure recording procedure required manual resetting of 
the recording pen and restarting of the recording trace every 750 sec- 
onds, and the changing of the recording chart every 3000 or 3750 sec- 
onds. All pressure data was recorded at a plotter pen speed of 1/50 
inches per second. As the length of the recorder chart was 15 inches, 

@ single sweep lasted 750 seconds, Usually it was possible to record 
four sweeps on the same chart page by lowering the pen position two 
inches lower for every sweep. ; 

Although the speed of recorder pen movement in the X-direction 
ef the chart was nominally 1/50 inches per second, it was found to vary 
slightly between the two recorders, For a better data control, a low 
voltage source was connected through a switch to the input terminals of 
both recorders, Whenever the switch was closed, instantaneous displace- 
ment of the recorder pen from the normal pressure trace provided a 
common time marker on both charts, 

The procedural steps discussed in the various paragraphs were 


interspaced with the actual pressure data acquisition steps, Excluding 
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these steps the outline of the actuel pressure measurement procedure. 
can be summarized as follows: 

(1) Immediately prior to the stert of a test run, the transducers 
were connected through the manifold to monitor pressures at any two of 
the twelve sensor locations in the flume, 

(2) The measured differential pressurs between the flume and the 
reference source was observed on the transducer meter dials and was ad~ 

justed to zero by raising or lowering the water level in the reference 
| reservoir with a hypodermic needle. 
| (3) The clamps on the pressure line from the reference reservoir 
- to the manifold were removed. Al] pressure lines connecting the sensor 
abe 2 the flum he manifolds ware clamped tight. This exposed 
both chambers of the transducers to an identical pressure. 

(4) The experimental run was started by activating the test 
chronometer and by the recording thse zero differential pressure output 
from the pressure transducers, A record of the zero differential pressure 
was essential as a common reference in data reduction and for correlation 
of pressures measured by the two independent transducer systems, A cone 
tinuous output of 30 to 60 seconds in duration was recorded, 

(5) Recorder input was interrupted by turning off the protector 
switch, 

(6) Preeselected sensor pressures line clamps were rsmoved and 
replaced on lines connecting the manifold to the reference source reservoir. 

(7) After a delay of approximately 20 seconds to allow manifold 
and the transducer chamber pressure to equalize with the pressures at 
sensor positions in the flume, the protector switches were turned on, and 
simultaneously the chronometer (time) reading was noted and voice recor- 


ded on tape, The identifying numbers of the sensors sampled were elso 
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voice recorded at this time. 

(8) When the short period pressure sampling procedure — 
employed, a pressurerecord of 50 to 60 seconds in duration was recor- 
ded. Otherwise, a continuous recording of pressure from a Single sen- 
sor was maintained. 

(9) At the end of the sampling period the recorder input pro- 
tector switch was turned off, 

(10) Clamps from two other sensor pressure lines were removed 
and replaced on the lines from the sensors already sampled, 
(11) Steps 7 to 10 were repeated until the end of the experiment. 

555 Integrated Data Acquisition Procedures 

The different procedures described previously for bed-wave 
-geometry, velocity, and pressure measurements were integrated together 
in a typical experimental run. Generally, the depth-wire readings and 
velocity measurements were carried out during the recording periods for 
the pressure data, Various other observations, bed-wave shape mainte. 
nance in the observation area, placement of time markers, etc, were also 
made during these periods. 

All data acquisition tasks, including visual observations and 
manipulation of equipment were performed by the principal observer who 
was positioned at the observation window. All equipment was located 
within reach from this position. The assistant observer, whose duties 
were to restore two-dimensional bed-forms at the screens and to maintain 
the shape of the bed-forms in the approach channel, was positioned 
approximately six feet upstream of the observation area. Both observers 
maintained their positions with minimum movement throughout a test run 


to avoid disturbing the flow in the flume, 
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Considerable experience on tho part of the cbserver in integers. 
ting all of the procedural steps and other tasks was required to success= 
fully obtain experimental data, The amount of data obtained was also, 
to a large degree, restricted by the numerous tasks that the observer 
had te perform in the available time, 

All observations were tied to the time base of each test run. 

At the beginning of an experiment, a digital chronometer reading in sec- 
onds, was activated and all times of all subsequent observations of bed 
“wave geometry, hydrogen bubble photographs and pressures were noted and 
recorded, These times provided the common reference through which ine . 


stantaneous positions of various probes in the flume and consequently 


geometry could be determined. 

In addition to the boundary geometry, velocity of flow and 
pressure data, readings of water temperature, flow discharge and the 
absolute msan depth of water in the flume were recorded before and efter 


an experimental run. 


5.6 Data Reduction 

After the conclusion of a test run all data was reduced into 
digital form. 

All observations and data readings contained on the audio tape 
were transcribed first. Pressure sampling times and sensor identifiers 
were marked directly on the pressure record charts, Bed-wave profile 
depth-wire readings and tho mean time of scan wers recorded in tables, 
During the initial data transcription only the times, wire data and film 
frame numbers of the hydrogen bubble photographs taken during the experi- 


ment were transcribed. 
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Depth-wire cbservetion data was transcribed directly from the 
audio tape into two tables, In the first, the mean time of the scan (from 
the start of the experiment) and individual wire readings were recorded. 
In the second table, the measured or estimated positions and elevations 
of the crests and toes of observed individual bed-waves were entered, 

The films of the hydrogen-bubble photographs were developed and 
prints 3.25 by 4.75 inches in size were made commercially. A typical 
hydrogen bubble photograph is given in Figure 5.1. Depending upon the 
size of the hydrogen bubble images on the photograph the digitizing of 
the visual data was made either directly from the film or from the photo. 
grephs, It was found that if the images of the hydrogen bubble rows 
were large the distances botween adjacent rows of bubbles could be 


scaled sufficiently accurately from the photographs. If, however, the 





Figure 5.1. Typical Hydrogen-Bubble Photograph 


Be 


ott: aan ehimenth, badhzeqana? Age sab yldayeendy, ertndiage 


Srrgt } eas 


noc fieve 


(2 4 .vileloirseeee ohn evew ante ab eatond £348 aad 2 | 


4 . - - . ol oneal uy f " > > J , odd 
es 3 on i = .& Vo» tf son re 68 anal ajab Laue “Bi; 


evened . FT e8tiy tndosodgq ei? worl qlotaiwons ctnotostive bela 















fae 7 : 
a 
Wt te omtt wane add ,fouk? et al otedind ont sind gat 98 
— evie Saubivitent dae. \temetconpee add iyo 
- z eraltiney hatentzae wo Semuswn, odd soiges teroogn add. et 
+e en ii ieiiacaliae taubbetbas heryeude 0 neat one aseeno aa? Te | 
sob wae sdasyyotertq. ofddud~aspegingt adi. to alt? off | 


— .{.2 eugtt af aeets,. ef dqatgotodg, elited geggamas 
hotahh ats seat rots eos o §enanl efddird aegoriyd oid To qane 


olddad senetiad edt to senant edd i dartt Bxtvol eaw JT ; alias 
bivoo eeldded to ewor Invoa ba neewded eoonnsetb af? egeal ‘ 


7 i fe ik 





87 


photograph contained rows of bubbles from 211 three wires, and the 
images were therofore small, the scaling from the photographs became less 
accurate because of the small absolute distances between the rows of 
bubbles and the digitizing was made using @ microfilm reader to ree 
the enlarged film image on 0,1 inch square grid, 

To digitize the hydrogen bubble information, the scale of the 
photograph, or of ths projected image, was determined first from the mean 
of the distances between the images of two adjacent hydrogen bubble wires 
scaled at two different levels below the surface of the flow. As the 
true distance between these wires in the flume was known, the scale of 


photograph or the projected image could easily be calculated. In general 


differences, up to one per cent of the total distance, were found to 
exist eon some photographs. These errors were attributed to a possibie 
stretching of the film or of the print, or possibly, to the combined 
optical effects of air, glass and water if the camera lens was slightly 
tilted from the perpendicular to the glass wall when the exposure was 
made. Averaging of the measurements, therefore, tended to reduce the 
error of the true scale of the photograph. 

Information to determine the instantaneous velocity profiles 
above the sediment bed, the location of a velocity profile in relation 
to a wire, and the flow depth at the wire was extracted from the hydrogen 
bubble photographs. The horizontal distances between two or more rows of 
hydrogen bubbles were accurately scaled at different levels between the 
sediment bed and the free surface, The number of levels at which the 
spacing between the rows of bubbles wes scaled varied from four to ten 


depending on the apparent uniformity of the velocity profile. The dis- 
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tance between the wire snd the upstream beunding row of bubbles for the 
determination of the space interval between the rows, was also measured. 
In addition, the distance from the free surface to the sediment bed at 
the wire was also scaled, | 

The transcribed information was summarized for each instanta- 
neous observation of velocity. This summary contained film frame number, 
time of exposure on the test run time clock, wire pulsation frequency, 
wire number, data to determine scale of photograph, depth measurement at 
wire, scaied distances from the wire to the bounding row of bubbles, 
sealed distances between rows of bubbles and the number of row intervals, 


and the scaled distances from the free surface to each level at which 


In the reduction of pressure data, the first step was the tran- 
seription from the audio tape of all anslog trace identification infor- 
mation of their starting times, and special marker times into the pres- 
sure data record charts, 

Next, a reference pressure common to all sensor prebes was eSe 
tablished on the original pressure record charts, For this the sensor 
calibration information obtained in the pre-test instrument check with 
no flow in the flume was used, The differences in the transducer out- 
put voltages between the applied zero differential pressure (both 
chambers of a transducer connected to the reference reservoir) and the 
constant hydrostatic head in the flume were determined for each sensor 
probe, The cormon reference pressure Cn the data charts was then 
established by indicating these differences above er below the zero 
differential pressure trace for each corresponding sensor. It must be 


noted that the zero differential pressure reference was used only as & 
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control to maintain true output relationships between the different sen- 
sor probes. For this reason, the level of water in the reference reo 
Servoir had to remain constant during the pre-test instrument check and 
throughout the experimental run. It could, however, be raised or 
lowared between these steps. 

All analog pressure data records exhibited amplitude fluctuations 
ranging from 0.001 to 0.005 inches of water, As the frequency response 
of the pressure measurement system precluded meaningful interpretation of 
the fluctuation pressures, only the mean pressure information was ex- 
ceaoted from the records, Uppsr and lower envelope lines were drawn to 
fheriae all but the extreme peak pressures recurring at intervals in 


excass of 30 saconds, The mean pressure line for each samp] 


was then drawn bisecting the distance between the enveloping lines. 

The actual digitizing of date was made by measuring the voltage 
difference between the appropriate sensor reference pressure line and 
the mean pressure. In converting the analog pressure data into digital 
form two criteria were considered, If tho mean pressure line was 
constant, only a single pressure value was determined at the mid-interval 
time of the sampling period. If, however, pressure variation was noted, 
pressures wore determined at the beginning and the end times of the 
sampling period and between, at intervals depending upon the pressure 
gradients but not exceeding 50 seconds. 

The pressure data were grouped for each Sensor probs separately, 
in tebles which identified the probe, its position in relation to the 
flume (i.e. its horizontal and vertical coordinates) and the times the 


voltage values of each observation. 
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CHAPTER VI 
EXPERIMENTAL RESULTS 


6.1 Introduction 

| The results of the experimental measurements are presented in 
this chapter, First, various pertinent terms and reference systems are 
defined then the descriptions of raw data treatment procedures are given. 
Criteria for the auality of data are discussed and data sets for final 
analyses are idontified. Final data treatment procedures are described 
and summaries of results are presented. Included in this chépter are 
also the results cf qualitative observations relevant to the experimental 
gGa1LS, , 

Exclusive of the preliminary trials a total of 29 individual 
experiments, under various flow and depth cenditions were performed in 
this investigation. 

Because of bed-form deformation, four tests were terminated before 
completion, and in the remaining 25 experiments, observations and measure~- 
ments of various parameters were made on 55 individual bed-forms. Because 
data contained more Pray 20,000 individual bits of digital information, 
all preliminary data treatment and plotting were performed by an IBM 


360/67 computer and Calcomp plotter. 


6.2 Bedeform Coordinate System _and Nomenclature 
The raw experinental data, transcribed from the audio tape and 


the analog charts, provided information about the instantansous absolute 
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values of the measured perimeters at the fixed locations of each probe 
and at the instant of time of cech observation, 

To establish the instantaneous positions of the various probes, 
and consequently the spatial variation of the measured parameters in re= 
lation to the geometry of cach bed-form, two rectangular coordinate 
systensem.ons stationary and one moving=-were employed. 

The origin of the stationary coordinate system was for conventence 
placed at the intersection of the first (upstream) depthewire and the top 
surfece of the instrument plate, The ordinate was aligned vertically 
along the wire, 4nd the positive direction of the abscissa was pointed 
downstream @long the flume, Al] observations were made in the stationary 
reference frame, 

The origin of a moving reference frame was placed at the inter. 
section of the vertical lino passing through the crest of a moving bed- 
form and the horizontal lines at the mean elevation (in the fixed re. 
ference system) of the toe of the bed-form upstream of the reference crest. 

The location of the origin ef the moving coordinate system for 
each bedeform was selected because crests and toes were the only clearly 
recognizable features of & bedeform. Tho vertical position of the moving 
origin at the toc, upstream of the reference crest, was chosen because 
the elevation difference between the toe and the crest downstream served 
&s & convenient reference parameter for non-dimensionalizing of the bed- 
form profile geometry data, A crest, by definition, was the intersection 
of the leo slope of a bed-form with the surface of its upstream face, 
Similarly, 2 toe was defined by the intersection of the lee slope with 
the upstream face of the adjacent downstream bed-fornm,. 


The coordinate systems and various descriptive terms are 
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illustrated in Figure 6.1. 

Bedeforms in @ train were numbered in their sequence of arrival 
at the observation area, The limits of a bed-wava were considered to 
extend upstream from the toe of wave (a) to the toe of wave (a +1). 

However, in evaluating bed-form wave lengths a greater variabil. 
ity was found between the wave lengths determined from toe to toe than 
from erest to crest measurements, This was caused by the cyclic sliding 
of particles on the lee side of a bed-form which advanced each tos in 
sudden steps cf random length and at random time intervals, To reduce 
data scatter, psrticularly in computing parameters that involvod the wave 
length of a bed-form, the adjacent crest to crest distances were adopted 
as the wave length measures, To differentiate from the bed-form defi- 
nition, the space between two adjacent crosts was called a range, A 
range wes identified by a number corresponding to ths number of the up. 


stream bed-fornm. 


6.3 Flow-Roundary Stability 





6.3.1 Prolimninary Observations 

The objectives of the expsriments were to obtain measured data 
from which the characteristics of various flow and sediment transport 
parameters within the bounds of a two-dimensional bed-form could be re. 
constructed. To obtain meaningful measurements, the characteristics of 
the structure of the flow and the response of the sediment particles had 
to rem&in essentially unchanged in respect to the geometry of the bed. 
form during the observation period, when the bed-form traversed across 
the area containing the stationary instrument probes, 

Although all precautions were taken to insure that bed-forms 


were stable during the observation period, some tests were terminated 
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before completion when the bed-forms under observation changed their 
heights or lengths, or in some other way became distorted from the re- 
quired two-dimensional form, This observation suggested that signifi- 
cant variations in the geometry of bed-form under observation may have 
been visually undetected in other experiments. As temporal variations 

of bed-form geometry during the experimental runs could not be monitored, 
a method was devised to evaluate individual bed-form stablity characteris- 
tics from experimental bed-form geometry data, 

6.3.2 Criteria for Bed-Form Stability 

To determine if a bed-form was in equilibrium during the period 
of its observation, and to insure that the related pei eeriey data were 
suitable for further analysis, 4 common set of standards was adopted to 
all bed-~forms in the experiments. 

The standards were established on the principle that a bed-form 
was in equilibrium if the flow and the sediment response characteristics 
in the space bounded by the upstream and downstream limits of a bed-form 
remained invariant throughout the period of observation. 

It was assumed that the hydrodynamic characteristics of flow were 
constant during the period of observation if the following parameters re- 
mained constant: (1) the mean velocity of flow entering and leaving the 
spatial element above a bed-form, (2) the wave length of the pedetorse 
and (3) the height of the upstream bed-form. 

Similarly, it was assumed that the sediment response characteris- 
tics were constant if the geometry of the bed-wave, as measured by its 
wave length and height, remained constant during the observation period. 

As the imposed flow discharge was constant throughout an experi- 


ment, @ test for the constancy of elevation of a crest above a reference 
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provided the test for constancy of flow velocity. The information re- 
quired for this test was available directly from the observed data, 

| The information required to test the constancy of bed-form 
heights w8s computed from the instantaneous observed elevations of 
toes and crests, 

Finally, the information required to test the constancy of bed- 
form wave lengths was obtained from the observed data of the instanta- 
neous positions of the two crests in the stationary reference frame in 
the flume. 

6.3.3 Statistical Analysis of Bed-Form Equilibrium 

Initially, an attempt was made to evaluate bed-form equilibrium 
by plotting crest elevations, wave heights, and crest positions on a 
time base, Because of the limitations of data acquisition techniques 
and natural effects in the flume, these plots showed considerable data 
scatter. To avoid inconsistencies in evaluating data by visual means, 
linear regression techniques were adopted to test stability of experi- 
mental bed-forms. In the regression analyses, crest elevations, wave 
heights, and wave lengths were considered as the dependent variables, 
The time of each observation was treated as the independent variable, be- 
cause its measurement incurred less error than measurement of other vari- 
ables, All computations, including the determination of regression 
coefficients, were performed by an IBM 360/67 computer using an avail- 
able library program. 

To test if a crest or a toe remained at a constant elevation 
during an observation period the regression coefficient, representing 
the slope of the regression line with time as the independent variable, 


had to equal zero. It was tested at a given significance level under 
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the null hypothesis H, : @= 0, against the alternative H : B £0 
by the statistic 


where b is the regression coefficient, and V(b] is the variance of its 
axpeeted value which is distributed as a Student's t.variable with 

(n - 2) degrees of freedom and where n represents the number of obser. 
vations (Keeping, 1962). 

To determine if the wave length of a bed~form was constant during 
an observation period, a joint statistical anslysis of the regression 
coefficients by and by of the horizontal positions of two adjacent crests 
with time as independent variable was carried out, The theoretivel back. 
ground of the statistical test which was used is given by Brownlee (1960). 
The null hypothesis under a given level of significance was He :By =Bo, 
against the alternative Hy +64 # Bo. The statistic for testing the nuli 
hypothesis was 


by - bo 


which is distributed as (ny 2) & (ny - 2) degrees of freedom where 


(6.2) 


Ny and Ny represont tho number of observations in data for regression 
lines 1 and 2 respectively. 

The selection of individual bed-form measurement data for further 
analysis was based on the results of statistical evaluation of the con- 
stancy of five bed-form and flow parameters during the observation period, 
namely: (1) submergence of the upstream bed-forn erest, (2) submergence of 


the bed-form crest, (3) height of the upstream form, (4) height of the 
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bed-form, and (5) wave length of the bed-form, Significance of all 
parameters was considered and data were accepted brfeondl Clonal if 011 
five parameters were non-significant at the 0.05 level. Otherwise, bed- 
form test data were accepted if three or less parameters were not 
significant at .02 level, or two or less parameters were not significant 
at 0.01 level, or not more than one parameter at 0.005 level, As a result 
of the abeve conditions only 17 out of 55 bed-form profiles were accepted 
for further analysis, A resume of the statistical bed-form equilibrium 


test results is given in Appendix C, 






6.4% Bed-Forn Profile Gaometr 

Stable bed-forms transleting downstream with constant velocities 
reflect the relative equilibrium between the hydrodynamic force and sedi- 
mont response systems, The geometric characteristics of their profiles 
may therefore serve in explaining some of the physical aspects of the 
entrainment. and transport of sedinent as bed~load along bed-form surfaces, 

A summary of the basic bed-form and flow parameters of bed-forms 
selected for final anslysis is presented in Table 6.1. 

The mean depth of flow ad above a bed«form was determined by esti. 
mating the area occupied by water between two adjacent crests and sare 
it by the mean wave length of ths bed-form, The flow area was estimated 
by the trapezoidal rule methed dividing the wave length into ten equal 
parts, 

The mean depth of flow above a crest, daa or doo were determined 
from the difference between the elevation of the free flow surface and the 
mean elevation of the crest measured in the stationary reference frame. 

The mean height of & bed~forn Hy or Hp were determined by sub- 


tracting the mean elevation of the tee from the mean elevation of the crest 
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of the bed-form 2lso measured in the stationary reference frame, Ths 
mean reference height A, of a bed-form was calculated by subtracting the 
mean elevation of the origin of the moving reference frame from the mean 
elevation of the reference crest, measured in the stationary reference 
frame, 

The mean wave length L and the mean velocity of translation Vp of 
a bed-form were determined by the joint regression technique described in 
Section 6.3.3. 

The representative profile geometry of an individual bed.form was 
determined in three steps, First, the bed profile coordinates originally 


observed in the stationary reference frame were transferred into the moving 


Second, sll of the observed surfece coordinate points were nonedimension- 
alized and, finally, a mean profile for a bed-form was calculated from the 
dimensionless ccordinate data. 

To compute bed-form surface coordinates in a moving reference 
frame the coordinates of the origin of the moving reference frame xp and 
ho for every observation time had to be established. The vertical coordi. 
nate h, for sach bed-form range was determined from the mean elevation of 
the tos contained in the appropriate bed-form range. The horizontal 
coordinate x, was determined aither from observed data or was estimated 
from regression coefficients relating time and the position of the reference 
erest in the flume. The application of the regression coefficients was 
found to be useful in extending the utility of some data points when the 
governing crest moved outside the observation area, 


The horizontal and vertical bed-form surfacs coordinates X, and Hs 
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in the moving reference frame were computed and nonedimonsionalized in 


accordance with the following equations: 


Xs = ag. Xo (6,3) 

Hy = hs = hy (6,1) 
x 

x= Xs 6. 
= (6.5) 
H 

H = — (6.6) 
He 


where x, and h, represent the instantaneous horizontal and verticaal 
coordinate of sediment surface data points in the staticnary reference 
frase. 

Finally, the representative profile of an individual bed-form 
was determined by grouping 811 available dimensionless bed surface data 
points into nine segments along the bed-form length, 0.1 L long each, 
and computing their mesn position in coordinates X and H in each segment 
in the dimensionless coordinate system according to 
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(6.8) 


where mis the number of coordinate points in & segment. 

Tho representative mean profiles of bed-forms selected for analysis 
are presented in Figures 6.2(a) to 6.2(a), Multiplying factors L and He 
which may be used to convert a dimensionless profile into natural coordi- 
nate system and the relative size of particles Deo are also shown for 


each bed-form, A bast fit smooth line was drawn between the weighted 
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Figure 6.2(a), Mean Dimensionless Bed~Form Profiles, 
Test ll, Range 2, and Test 12, Range 2. 
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Figure 6,2(c). Mean Dimensionless 3ed-Form Profiles, 
: Test 20, Range 2, and Test 20, Range 3. 
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Figure 6.2(d). Mean Dimensionless Bed-Form Profiles, 
Test 21, Range 2, and Test 2], Range 3. 
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Figure 6,.2(6). Mean Dimensionless Bed-Form Profiles, 
Test 22, Range 2, 
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Figure 6,.2(¢). Mean Dimensionless Bed-Form Profiles, 
Test 26, Range 2, and Test 26, Range 3. 
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surface coordinate data points. The relative size of sediment grain 
shows that in the vertical the surface data points generally fall within 
ene half grain diameter from any estimated surface profile line. This is 
within the predicted dopth-wire resding error. 

A Soden bed-form profile, developed by averaging H values at 
each 0.1 X, is presented in Figure 6.3. Also shown are the limits of 
scatter of H values at cach X. 

The average horizontal position of the bed-form toes was found 
to be at X = -0.935. With the exception of one bed-form, the scour below 
the ae toe elevation in the trough region was less than one particle 
diameter deep. 

The smooth profile of the back of the bed-forms had its maximum 
slope in the range -0.60 < X < -0.45. However, on five separate bed~ 
forms (12/2, 17/2, 17/3, 24/2, and 24/4) the smooth profile appears to 
terminate in the first 0.1 L segment upstream of the crest, followed by a 
sharp rise to the crest. This effect was observed visually during sone 
of the experiments when it was noted that often a single particle at the 
apex of the crest rested slightly above its neighbouring particles 


upstream. 


6.5 Entrainnent and Transportation 

6.5.1 Genoral Observations 

Visual observations indicated that the movement of individual 
bed particles as bed-load over bed-forms was intermittent, Particles 
were entrained on the upstream side of the bed-forms, transported along 
the surface in steps of random lengths and deposited in the lee of the bed- 
form crests. The particles then remained at rest buried within the bed- 


form until the continuous entrainment and transport of particles from the 
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Figure 6,3. Average Dimensionless 3ed-Form Profile 
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upstream side translated the entire bed-form downstream and re-exposed 
the buried particles for a repeated entrainment, transportation, and 
deposition sequence, 

~The poriod of time during which a particle remained at rest when 
buried in ths bed-form depended upon the position at which the particle 
was depos ted on the lee slope downstream of the crest. The further 
below the crest it bectme buried the longer was its time of rest. The 
distance of its path to the crest after re-exposure and entrainment 
was similarly longer if the particle was buried near the toe than near 
the crest. 

Viewsd from a stationary reference frame an individual particle 
is transported dewnstream along an irregular sawtooth path shown in 
Figure 6.4, However, when viewed from a reference frame moving at the 
speed of the translating bed-form, the paths of tho particles describe a 
near triangle, All particles buried within the bed-form recede in the 
upstream direction at a constant velocity of bed-form translation, When 
exposed on the upstream side of the bed-form they catch up with ths crest 
with an average velocity much higher than the bed-form velocity itself, 
If the identity of individual particles is ignored and the bed-forn 
material 4s considered as a deformable continuum, its motion may be re- 
presented by a system of "streamlines" illustrated in Figure 6.5. It 
suggests that the material in the bed-form rotates about a center located 
at the crest of the bed-form. Because of the translation of the bed-forn, 
a relatively high angular velocity of rotation of material persists near 
the crest, The observed paths of particles partially confirm this 
phenomsnon. Particles arriving along the bed were observed to lift 


slightly above the bed immediately upstream of the crest, and, after 
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Figure 6.4, Particle Path in Stationary Reference Frame 





Figure 6.5, Particle Path in Moving Reference Frame 
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rounding the crest, were deposited nosr the crest on the les slopes of 
the bed-form. 

6.5.2 Distribution and Rates of Entrainment 

Particle entrainment rates (scour) and distribution and, hence, 
the production of sediment for transport along bed-form surfaces wers 
estimated by two methods, namely: (1) depth-wire method, and (2) profile 
translation method, 


In the depth-wire mothod, local rates of particle entrainment 
were calculated from all available depth-wire readings of a bed-form 
surface profile from the changes cf bed elevation at a depth-wire between 


two suksequent observations using the following formula: 


hi 





where hy and hy are the depth-wire readings at observetions times ty and 
ts respectively. The horizontal position of the computed entrainment rate 
wes considered to be st mid-point between the two positions of the depth- 
wire in the moving reference frame as iliustrated in Figura Goce 

The Brateivation of entrainment rates was obtained by grouping the 
computed local entrainment rate data points into appropriate 0,1 ‘8 
segments along the bed-form profile and estimating the mean rate of 


entrainnent from tho following formula: 


Nn. : 
Di: 
yes Se (6.10) 
Tig 


where ng is the number of data points in each segment. 
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Figure 6.6. Definition Diegram for Entrainment Rate 


The volumetric rate of bed-material entrainment g, per unit 
width of bed-form from each segment was calculated by the following 


formula: 


g = O1Lh (6.11) 


Tho distributions of typical entrainment rates h', calculated 
from results of tests 28/2 and 28/3 by the above mathed and rendered 
dimensionless by multiplying by the ratio of bed-form period to bed-form 
height, are illustrated in Figures 6.7(a) and 6.7(b). The horizontal 
position coordinzte X has been non-dimensionalized in accordance with 


equation 6.5. The scatter of the data reflects both the error in measure- 
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Typical Instantansous Entrainment Rate Results 


Figure 6,7(b). 
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ment and the natural fluctuations of the sediment surface, due to inter- 
mittent movement of particles at the same relative position on the bed- 
form as it translates downstream, The average rates of entrainment h' 
distribution along the bed-forms are also shown. 

In the second method, the volumetric rates of particle entrainment. 
and distribution from the surface of a@ bed-form were determined employing 
the mean bed-form profile data and the velocity of its migration down- 


stream Vp from 
Eb = Vp(hy oo h5) (6.12) 


where hy and h” represent the mean elevations of bed-form profile at the 


* 
2 
upstream and downstream bounds of a segment 0.1 L long as illustrated in 
Figure 6.8, 

The experimental entrainment rate and distribution results for 
the seventoon bed-forms are presented in Figures 6.9(a) to (q). The re- 
sults of both methods of calculation are given for comparison purposes, 

Entrainment of particles occurs at all locations of a bed-form 
where its surface slope is greater than zero. Scouring processes are 
therefore active downstream of the deepest position of the trough and dé= 
position occurs between the trough and the toe upstream. 

The maximum entrainnent rate, as estimated by the moan profile 
method, occurs when -0.6 ¢ X < -0.45 which corresponds to the location of 
the maximum slope of bed-form surface, From the maximua position the rate 
of entrainment drops toward the crest and, in the opposite direction, to- 
ward the trough of the bed-form. In general, the same result is obtained 
by the depth-wire method of calculation but because of the numerical 


averaging procedure with a limited number of data points in each segment, 


ie 









a 


a ¢ * * & a -— eas Zi er as 


} it oi 3 bets Suan 


fe 38 eae 


‘ ‘ f j ee . 2 mets “is saotss 


; 4 = ° 
ae t Me i 
: “i@ } f¢ 
MP Fae 
j 
: 
A 
: % 4 ; 
- > 7 
» 
, 
ie 4 
- id 
5 Gc 1 . 
- ° 
f 
’ _ 
“ lentxG 
* 
. i ; - = 
ng = - . - ; ’ a 
ne . < nid £7] y a 
i ‘ 
* a - . i _ a 
t tue 7 ry ‘ MmoVwS do ENIFOOG 
* 
wae or > is, f 7 4 ; ror: 7 at . J a _— : 
et Pea seal + MUSSELS SS ES ,OFes sean le ine susataan @fT- 
3 . JL é 
3 7 


— c = "ae, sy PRS af Pai B ba a as ro o noLdvw P+ p= > > 4 - s.2 20~. wert ad W700 
ats mt elt. aolkdts < tps ite nm ~~ . if 
' 4: M 4% gh. O'S anak F _seoepem cee abod, Xe he a ts 


utabtsogth otteoa att et pews odd-brayot eqotb does 
fase | 


= 

7 _ ee oe 

wy : Tas tT | wer. ( . ‘ As - 
—s a . . t 


et 
an 
_ biaehahetentiads 









118 


the entrainment distribution is not as uniform as the distribution by the 


profile mothod. 





Figure 6.8. Definition Diagram for Entrainment Rate Calculation by 
Profile Translation Method 
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The distribution of sediment flux per unit width of a bed-form 


was then calculated by 


CLs = > (gy )5 (6,15) 


where Gus is the volumetric rate of sediment flux at the downstream limit 
of segment s as measured by the volume of bed-fora scoured, 

To determine the distribution of sediment flux on the surface of 
a bedeform from particle entrainment rates required also the knowledge of 
the position of zero net transport. In the present calculations this 
location was assumed to be at the moan position of the maxinmm pressure as 
determined from pressure measurements, and was verified by the total rate 


of sediment flux per unit width of the bed at crest G,, computed from 


Semen ria eo a (6.25) 
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Figure 6.9(a). Distribution of Particle Entrainment and 
Transport Rates, Test 11, Range 2, 
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Figure 6.9(b). Distribution of Particle Entrainment and 
Transport Kates, Test 12 Ranges 2, 
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Figure 6.9(c). Distribution of Particle Entrainment and 
Transport Rates, Test 17, Range 2. 
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Figure 6.9(d). Distribution of Particle Entrainment and 
Transport Rates, Test 17, Range 3, 
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Figure 6.9(g). Distribution of Particle Entrainment and 
Transport Rates, Test 21, Range 2. 
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Figure 6.9{h). Distribution of Particle Entrainment and 


Transport Rates, Test 21, Range 3. 
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Figure 6,9(1). Distribution of Particle Entrainment and 
Trensport Rates, Test 22, Range 2 
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Figure 6.9(n). Distribution of Particle Entrainment and 
Transport Rates, Test 28, Range 2. 


nas q 
Aube mr Deke 


—_ 





ms 









A im 0S 
Depth wire method, gp, 


Profile translation method, G. 
Profile translation method, ¢, 


Goo 


@®pPr do Bs 


Gog x 10% (cu in/see)/in 


a6) 
i Oho Ue Sauime fim Die Om Do Simm g lh maO yd mreO se a0 se 050 


x 


Figure 6.9(0). Distribution of Particle Entrainment and_ 
Transport Rates, Test 28, Range 3. 


OS i i ea | eet T i 
Depth wire method, G_ 


133 


(cu in/see)/in 


. alt 
&, x LO 










. a OE MUI PR ee 


| ei” hadion otty diged se 
bp ,hotfenr etiv fideo 

not! «banter cotta Tuners a Litort « 
| my. pattie netic fonds? alia 2 
a 


: of” 


Depth wires method, Gu 
, a o 
Depth wire method, ¢ 
Profile translation method, G,, 
Profile translation method, g, 


@DpeO & 


Ghe 
20 
18 
16 
Sc 
ee 
an tie 
& 
ce 
ee. . 
’ 12 
0 
9 10 ~ 
F ees 
pire CAR 
: 6 SPS 
i Aa 
y/ / 
4 fA e/ 





LTO P2059 076" 0.7 ~0.6 =-0.5 =O.4% -0.3 -0.2 -0.1 0.0 


x 


j ; Entrainment and 
Figure 6.9(p). Distribution of Particle 
e Transport Rates, Test 28, Range 4. 






& 


134 


&, x id (cu in/sec)/in 


ner 







Pp; en ee oe 


i gbesidot athe -diqed 
dt «GAT anit dget 
4 ot .( nen -netialenwrn? elect 

9 .b@len seftelsaer? olPiost 


od 


P| 


ec ros 





135 


ny hb <— S at aah hy, 3 2 oe m 
The sediment flux distribution on individual bod.forms is shown 


also in Figures 6,.9(a) to (q). 
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6.6 Hydrogen Bubble Velocity Measuromont Data 

6.6.1 Instantancous Velocity Profiles 

The term "instantaneous" as used in this context denotes a point 
velocity or 4 velocity profile determined from the separation distance of 
two rows of hydrogen bubbles which were separately introduced into the 
flow at a short, usually 0,1 second, time interval. 

The instantaneous velocity profiles parallel to the surface of 
flow were established from point velocities measured at different depths 
4n the flow above the bed, 

The number and spacing of depths at which point velocities were 


determined varied depending upon the uniformity of the separation distence 


0 


of the two rows of bubbles on a photograph. If the spacing betweon the 
rows was uniform bubble separation distances were determined only at the 
extrens limits of the uniform separation zone. However, where velocity 
gradients were pronounced, such as near the sediment boundary, separation 
distances between the rows were determined at additional depths to provide 
sufficient number of points to properly define the velocity profile. 
Although an attempt wes made to determine the spacing of the 
hydrogen bubbles at the top of particles, the reflection of light from the 
white coloured sediment obscured the rows of bubbles in the near vicinity 
of the boundary, Hence, the transcribed data from this region cannot be 
considered as accurate. In some photographs, however, because of more 
favorable viewing angle or shadows developed by slight lateral slopes of 
the bed, the bubble rows wore discerned at the top of boundary particles 


and the estimated flow velocities parallel to the bed were found to range 
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from zero to various finite values. The point velocity of flow at any 
depth below the surface was calculated by 
| (2, = a5)¥s 
Lvasine Tie cased” (6.15) 
Nye i 
where Uy 4s the estimated point velocity at a distance y above the bed, 
ay and a) ara the distances from the hydrogen bubble wire to the first 
(downstream) and second (upstream) row of bubbles respectively, n, is the 
number of full pulse cycles between the rows of bubbles used in computa- 
tion, T megeanents period ef one pulse cycle, andrgis the scale of thse 
photograph. 

: The instantaneous velocity profiles were determined from point 
velocities moasured at differsnt depths of the flow. The horizontal 
position of a profile was assumed to be at the mid-point between two rows 
of bubbles used in surface velocity calculation, Since the time of the 
observation was known, the horizontal coordinate of the moving reference 
frame X, was determined either by interpolating the position of the cone 
trolling crest between tio subsequent bed-form profile observations or 
from regression line of the crest position, and the coordinate of the ine 
stanteneous velocity profile in relation to the geometry of the bed-form, 


Nie was calculated in the dimensionless coordinate system by 
= — (6215) 
i 


where x, is the position of the velocity profile in the fixed reference 
frame. 


6.6.2 Flow Discharge Comparison 


The hydrogen bubble method of velocity measurement was vorified 
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against gross errors by comparing flow discharge measured by the flow 
meter with average discharge dcetemrined frox 1 all hydrogen bubble profiles 
in a test. The mean instantanoous flow velocity Up, was computed assuming 
linear velocity variation between adjacent point velocities on the 
instantaneous profile and using the local depth at the pulsing wire. The 
results of the comparison based on discharges per unit width of the flume 
are presented in Figure 6.10. Generally, the discrepancies between the 
observed and the metered discharges do not exceed ten per cont. They were 
attributed to two causes: (1) insufficient nurber of instantaneous velocity 
profiles in some tests to define accurately the average discharge, and 
(2) the retardation of flow near the side walls of the flume which re- 
duced the true discharge value along the centre of the flume. 

6.6.3 Dimensionless Velocity Profiles 

Preliminary oxaminetion of instantansous velocity profiles taken 
at short intervals of time and at the same relative positions to the bec- 
form showed, however, that large scale velocity fluctuations existed 
throughout the depth of flow. To develop average velocity profile 


- 


characteristics corresponding to the actual flow above an individual bed~ 
form by statistical means, 4 much larger number of instantaneous velocity 
profiles was therefore required than was avajlable. 

General considerations discussed previously in Section 4.3.4 
indicated, however, that the bed-forms were exclusively in the dune 
phase. Assuming that the flow characteristics above dunes were similar, 
average velocity profiles were developed by combining all available veloci- 
ty measurement information. 

ALL instantaneous velocity profiles were non-dimensionalized into 


two sets. In the first set, dimensionless velocities were. obtained by di- 


viding the point velocity at any depth of flow by the ins Se rdenaone mean 
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Figure 6,10. Comparison of Discharge Measurement. 
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velocity determined from the instantaneous profile. In the second set, 

the observed velocitics at different depths were divided by the instanta- 
neous velocity at the free surface of the flow. The dimensionless depth in 
both sets was calculated by dividing the vertical distance from the sedi- 
ment bed to the point velocity depth by the total depth of flow at the 
Wire. 

All instantanecus dimsnsionless velocity profiles were grouped 
by position of observation relative to the crest inte ten 0.1 L segments 
along a dimensionless bed-forn as previously was done with bed profile 
data, <A group of instantaneous dimensionless prefiles based on the in- 
stanteneous moan velocity and belonging to the segmont -0.3<X < -0.2 
4s illustrated in Figure 6.11. 

The mean dimensionless velocity profile in each segment was de- 
termined by evereging the instantarcous u/U, or u/Ug velocities at ten 
dimensionless depths, Starting at the surface of the sediment particles, 
the instantaneous profiles at 0.0625 4, 0,125 4 and every other 0,125 : 
level to the surface of the flow. 

The two sets of dimensionless velocity profiles of flow above the 
bed-forms are presented in Figures 6.12 and 6.13. 


6.6.4 Equal Moan Velocity Li 





The dimensionless mean velocity profiles given in Figure cole 
sllustrate local flow features at different locations above a bed-form. 
oh oe they cannot be used to demonstrate relative velocity characteris. 
tics elong tho length of a bed-form because each profile is referenced to 
the local depth of flow which varies along the bed-forn. 

The mean flow characteristics as represented by equal velocity 
lines along a bed-form are illustrated in Figure 6.14. The bed-form pro- 


file was doveloped from dimensionless mean profile coordinates given in 
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Figure 6,3 using L/H, of 10, and d.,/H, of 2.5 which were in the range 
of bed-forms investigated in the present experiments, The equal veloci- 
ties were calculated relative to the mean velocity at mean depth of flow 
above the average bed-form. 

6.6.5 Flow Turbulence Characteristics 

Examination of instantaneous velocity profiles clearly indicated 
that large scale turbulence fluctuations existed in the flow. Five in- 
stantaneous dimensionless velocity profiles all obtained from segment 
-0.3 < X< -0.2 in test 28 and illustrated in Figure 6.15 show, for 
example, that the velocities in the direction of flow at Y = 0.0625 varied 
from the mean by as much as sixty per cent. | 

To determine horizontal velocity turbulence distribution in the 
flow above the bed-forms from the available hydrogen bubble information it 
was assumed that characteristics of turbulence, like those of the mean velo- 
city of flow, were similar for flows above a given phase of bed-forms, 

Dimensionless turbulence intensities were determined from the 
standard deviation of the fluctuations of the dimensionless velocities 
u/U, about the mean in each segment of the bed-form and at the same dimen- 
sionless depth levels as were employed in the average velocity calculations. 
The dimensionless turbulence intensity results are presented.in Figure 6.16. 

For the reasons discussed previously in Section 6.6.4 the dimen- 
sionless turbulence intensities 0 presented in Figure 6.16 are indica. 
tive of local turbulence cheteotaristion only. The turbulence intensity 
variation in the flow above a bed-form is therefore illustrated in 
Figure 6.17. The dimensionless turbulence intensity is referenced to the 
mean velocity at mean depth UM and the bed-form and flow parameters 


L/H,, and = are 10 and 2.5 respectively. 
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6.7 Pressure Measurement 

6.7.1 Pressure Data Evaluation 

All pressure measurements were made et sensor positions stationary 
in the fixed reference system. To establish pressure distribution 
characteristics over individual moving bed-forms, the horizontal coordi~ 
nates of each sensor in the dimensionless moving reference system were come 


puted by equation (6.16) 
5 a Miee—eeress (6.16) 
L 


where iy, is ths position of & sensor in the stationary reference frame, and 
X_ represents the horizontal coordinate position of the origin of the mov- 
ing reference frame at the instant of a pressure observation, 

Typicel pressure moasurement results obtained by the continuous 
monitoring method end the intermittent sampling procedure are illustrated 
in Figures 6,18 and 6.19 respectively. The recorded pressures represent 
pressure differences msasured between the probe in the flume and the con- 
stant level reference source, Different pressure probes are repressnted 
by distinct symbols. 

The results given in Figure 6,18 show a general scatter of pressure 
data points in the order of 70.001 inches of water. Two points identified 
in the plot, however, bear special mention. Both points appear below the 
general trend of pressure distribution values at their respective locations 
in the bed, The time of observation of these two pressure points, however, 
was found to be the same indicating that the recorded drop in pressure of 
0,002 inches of water resulted from an external effect such as @ long 


period surge in the flume. 
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Because Similar long psriod pressure fluctuations were observed 


in pressure observation plots in othor tests, all selected test data wore 


reviewed to isolate the external effects upon the normal pressure varia- 


tions.. To identify externally caused pressure fluctuations, recorded 


‘) 


2 


pressure traces from two or more sensors were compared noting discrepancies 
from the normal trends at identical timos of observation and different 
positions relative to the bed-form. When an external pressure pattern 

was identified, its approximate amplitude and duration were estimated from 
two or more independent sensor output traces. A correction was then 
applied to restore the pressure record to the natural pattern. 

Typical results of pressure measurement by the intermittent sam- 
pling technique sre illustrated in Figure 6,19, Although considerable data 
scatter is indicated, a thorough analysis of the data disclosed that the 
scatter was due partially from oxternally produced pressure fluctuations, 
partially from small losses of water from the system, and possibly from 
small variations in the roughness of the bed in the epproach channel. 

Jo determine and, if possible, to climinate the effects of pres- 
sure data scatter dus to various external causes, and to establish pres- 
sure variation charactoristies due to flew above bed-forms only, a stand- 
ard prossure data evaluation procedure was evolved, The procedure was 
based on the fact that pressure at a fixed position in an equilibrium 
bedeform remains invariant, Hence, pressure variation along the bed-forn 
and at a fixed distance below somo reference elevation must retain indi- 
vidual and constant characteristics as the bed-form translates downstream, 

To evaluate the external pressure effects all sensor probes were 
classified into three groups according to low, medium or high elevations 
of their ports in the flume. AL] available pressure data from individual 
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sensors bslonging to each group were replotted and compared, Because of 
probe spacing along the flume, pressure observations at a point by differ. 
ent sensors wore time delayed and differences between the pressures roe 
corded at a point along the bed by different sensors were indicative of 
external factors influencing the pressure, Time Gependent tronds of 
external pressurs factors could then bo established separately in all 
groups and verified by comparison between groups, If these trends were 
consistent in time and magnitude, mean pressure distribution characterise 
tics were established by superposition of individuel pressure sensor out- 


put data in each group. 





. 


The results showed that, except for the ganeral scatter of data 
which could not be elimineted by timo trend analysis, the pressure diss 
tribution characteristics established st different levels within the 
sediment bed along a bed-form wore identical, Henee, within the precision 
of measurements the mean pressure distribution deta indicated no vertical 
pressure gradients in the bed, and all data from different lsvels were 
therefore employed to establish msan pressure distribution characteristics 
@long each bedeform as illustrated in Figure 6.20. 

The results of pressure measurement were non-dimensionalized in 


accordance with equation (6.17) 


(Gaz) 





whera Pol 4s the reference pressure at the origin of the moving reference 
frame, p is the mean observed pressure at any position in the bed, V 5 is 


the mean velocity of flow at the crost of the upstream bed-form, and g 
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Figure 6.21. Dimensionless Pressure Distribution, 
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Figure 6.21. (Continued ) 
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Figure 6.21. (Continued) 
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represents the gravitational acceleration, The dimensionless pressure 
distribution characteristics are presented in Figure 6.21, No pressure 
measurement results are reported for test 22 ranges 2 and 3, and test 
26 range 2 either because only partially complete data were obtained 


or external pressure effects could not be resolved, 






< ° . Pj2 i" if o é ) LAT t batt iL 
_ 
som em aolsadsatoar ans oak tod at 
~ ’ 2" — j : a atti 


a BB 


oes | Laer r0txe 10 





CHAPTER VII 
ANALYSIS AND DISCUSSION 


Fed Outiine 

The pietinental findings are analyzed with a specific objective 
to establish the relative contributions of different forces in detaching 
particles from the surface of active bed-forms. 

| First, the observations and measurement results are analyzed 

and the characteristics of bed-forms and flow, pertinent to the objectives 
the investigation, are identified from the experimental data, Then, 
force requirements to entrain particles from different locations along 
the surface of bed-forms are determined, and forces availeble for 
particle entrainment are estimated from the experimental results, Final- 
ly, relative contributions of flow shear and pressure gradient forces for 


particle entrainment are established. 


7.2 Analysis of Observations 
x AIOE AES SES ED 


7e2el Bod-Form Geomet 





Observations and profile measurements show that a typical bed-form 
may be divided into two regions distinguished by either predominant scour 
or deposition of surface particles, The region from which particles are 
Scoured extends from the deepest location of the trough downstream to the 
crest of the bed-form. It corresponds to the region ere the bed-forn 
surface curve slope is positive. The region in which deposition of parti- 


cles takes place extends downstream fron a crest to the deepest position 
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in the trough of the bed-form., Throughout this region the surface slope 
is negative, 

Each of the two regions may be further subdivided into zones of 
distinct bedeform surface geometry characteristics, In the region imme- 
diately downstream from the trough the bed-form surface profile is cone 
eave upward. However, the surface curvature becomes concave downward 
at a point of inflection, which from the various bed-form profiles was 
found to occur in the region -0.60<X<-0.45, In the majority of pro} 
files, the surface curve approached the crest smoothly. However, sur- 
face profiles of five bed-forms (12/2, 17/2, 17/3, 24/2, and 24/4) 
showed discontinuities in the segment upstream of the crests, This ree 
sult confirmed visual observations that a single particle often rested 
slightly above the adjacent particles at the apex of the crest. 

The rogion of particle deposition, downstream from the crest to 
the trough, was also found to have two distinct zones: {1) the lee side 
of the bed-form upon which 411 particles arriving from upstream were de- 
posited, and (2) the zone downstream from the toe to the deepest location 
in the trough. 

Observations showed that on the lee side all particles were de- 
posited near the crest. This inereased the slope of the lee side of the 
bed-form and when the critical intergranular shear stresses were exceeded 
the particles on the slope avalanched downward. Inmediately after the 
slope failure the lee side slope was less steep than the critical, But 
as the process of particle deposition near the crest was not interrupted, 
avalanching occurred again watth time the slope increased to some critical 
value, A statistical analysis of 324 individual lee slope angle obser- 


vations disclosed that the average slope of the lee side was 344°, rang- 
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_ ing from a minimum of 26° to a maximum of 45°, The average angle was 
higher than the angle of initial yield,>, = 31°, However, the angle of 
the slope after evalanching was found to be equal to the angle of repose, 
dy = 26°, Especially significant is the observation that the lee side 
could reach a maximum slope of 45°, A possible explanation is that 
there was inflow into the lee side of the bed-form increasing the effec. 
tive stresses between particles and helping to resist sliding at steeper 
Slopes éngles,. 

For deposition in the toe to trough zone, particles were supplied 
partielly by the transport in the upstream direction from the trough re- 
gion, and partially by some loose particles during the avalanching of the 
lee side of the bed-form. 

The distribution of the mean entrainment rates of particles from 
an equilibrium bed-form depends upon the local slepo of the bed surface. 
The entrainment rate was zero at the dsepsst position of the trough where 
the bed surface slope was zero, and was maximum at the inflection point on 
the upstream slope of the bed-form where the surface slope was maximum. 
The entrainment rate, therefore, increased from the trough to a position 
at approximately X = ~0.55, and then decreased toward the crest. 

The local flux of particles along the surface of the bed-forms was 
computed from the local entrainment results assuming zero net transport at 
the separation streamline reattachment position, The flux of particles, there- 
fore, increased in the downstream direction and was maximum at the crest, 

Since in these exporiments all particles arriving at a crest were 
trapped in its lee, the flux of particles moving past the crest had to equal 
to the rate of particle deposition in the lee of crest, The flux at the 


crest, determined from entrainment information, was therefore compared with 
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the particle transport estimated by the crest velocity and height. 
Generally the results were found to be in good agreement, and since the 
variations between the two results showed no consistent trends, the dis- 
erepancies were attributed to experimental nature, especially in estab- 
lishing the true position of zero net transport. 

However, the reattachment position in all bed-forms was located in 
a zone of finite positive surface slope and entrainment of particles 
occurred between the trough and the reattachment point. The total rate of 
particle entrainment was therefore in excess of the transport at the crest, 
as the particles upstream of the reattachnent point were not included in 
eos transport calculation. Small part of these particles were transported 
upstream and were deposited on the ped downstream from the toe. In the 
trough area, however, the bed appeared to be very loose, And it is 
possible that, because of the methed of measurement, the volume of particles 
entrained from this area, when measured in terms of volums occupied when 
deposited in the lee of a crest, was much less than predicted. 

7.202 Velocity end Turbulence 

The results of the mean velocity measurements have been presented 
in two forms: (1) by the mean local dimensionless velocity profiles Um and 
5. veS. Y given in Figures 6.12 and 6.13 respectively, and (2) by the equal 
velocity lines representation of a particular bed-form geometry--depth of 
flow condition given in Figure 6.14. 

Each of the mean dimensionless velocity profiles Un ves. Y 
summarize the local flow conditions at different positions relative to the 
geometry of a typical bed-form. The results show that the flow between 
two adjacent crests may be divided into regions upstream and downstream 


from the separation streamline reattachment position located in the segment 
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Bo.7 2X <~0.6. Upstream from the reattachment position, in the region 
Boros x < 0.7, a shear flow zone, characterized by large velocity gradients, 
was found to exist at depth levels 0.1<Y<0.5, Near the sediment bed, at 
depth levels 0.0<Y<0.1, the mean flow was found to be in the upstream 
direction toward the adjacent toe. In the flow above Y = 0.5 the dimen- 
sionless velocity gradients (aU,,/a¥) increased in the downstream direction 
and simultaneously extended toward the free surface of the flow. The 
maximum velocity gradients near the free surface in the flow occurred at 
X = -007- 

In the region downstream from the reattachment point the local 
velocity profiles became fuller in the downstream direction indicating 
development of boundary layer. The velocity gradients in the upper flow 
gones gradually decreased and the flow velocities became more uniform. 
Simultaneously the velocity eredients became more pronounced and were 
concentrated in gradually narrowing depth zones close to the sediment 
surface, The profiles showed that the mean local flow, represented by 
Un = 1,00 occurred in the depth zone 0.347404 at all positions along 
the bed-form except immediately upstream of the crest where the mean flow 
velocity occurred at Y= 0.24. The results also show that the velocity 
profile at this position filled out suddenly from X = -0.15 and that the 
flow above Y = 0.24 became nearly uniform, as the velocity at the free 
surface exceeded the mean local velocity by only four per cent. 

Equal velocity lines, based on the average flow velocity at mean 
depth of flow above a bed-form and shown in Figure 6,14 for a bed-form 
with H/L = 0.1 and do (fH,= 2.5, illustrate flow characteristics along the 
bed-form. It was found that the depth position of the mean velocity Uy 


varies considerably above the bed, being at Y = 0.1 immediately upstream 
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of the crest and at —— 0.4 and X = ~0,62, Furthermore, a cell of re- 
latively high flow velocity was found to exist immediately downstream of 
the crest. As a result, in the given example, the free surface flow 
velocities ranged from a maximum of approximately 1.35 iy at X = 0,95 

to a minimm of 1,12 Ung at X = -0.6, where Uy is the mean velocity at mean 
depth. 

The basic turbulence characteristics, represented by the standard 
deviation of the fluctuations of the dimensionless velocities u/U,, and 
summérized in Figure 6.16, show that in terms of local mean velocity Up 
the maximum intensity occurred in the trough region of the bed-form at 
X = -0.85 and at the depth level Y = 0,125 above the bed, and decreased 
toward the bed and toward the surface of the flow. Upstream from this 
position near the bed the turbulence intensity decreased rapidly, and in 
the vicinity of a toe at ¥ = 0.0625 the intensity was found to be 0,01 
4ndicating relatively negligible turbulence, Similarly, turbulence 
sntensities decreased in the downstream direction. However, from approxi- 
mately X = -0.60 the intensity at ¥ = 0.0625 was found to be larger than 
in the flow further from the bed. It decreased aware the central core 
of the flow and appeared to be minimum at Y = 0.50, but increased again 
toward the free surface. 

The distribution of turovlence intensities above a specific bed- 
form as represented by equal intensity lines are shown in Figure Oeics 
The maximum intensity occurred in the core recion of the roller downstream 
from a crest at approximately X = -0.85. Alow intensity region was 
found to exist throughout a central depth zone downstream from X = -0.70 


to beyond the crest of the bed-form. A zone of relatively high turbulence 


was found to exist throughout the full depth of flow above the trough 
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region which apparently limitedthe extent of the central, low-intensity 
turbulence region. In general, the turbulence intensity characteristics in 
the flow near the bed in these experiments agree qualitatively with the 
findings of Raudkivi (1966) on stabilized bed-forms, and Rifai and Smith 
(1971) on triangular elements. However, neither of these investigators 
found increasing turbulence intensities toward the free surface of the flow, 
or detected the zone of relatively high turbulence intensity penetrating 

to the surface from the trough region below. 

The results also show that downstream from X = -0.55 and near the 
bed, the distance of equal turbulence intensity lines oscillated from the 
boundary. Turbulence intensity appeared to be suppressed at X = -0.55, 
however it increased in the downstream direction and reached @ maximum at 
Xx = -0.35. The position of the turoulence suppression point corresponds 
to the lecation of the maximum pressure gradient and with the maximum 
positive surface slope of the bed=form. It is possible, therefore, that 
inflow of fluid into the bed occurred in this region which, according to 
a summary presented by Watters and Rao (1971), would have tended to decrease 
the turbulence level in the flow field above. Tho increase in turbulence 
intensity from X = -0.65 to X = ~0.35 appears to be due to turbulence 
generated in the developing boundary layer. 

The velocity and turbulence measurement results show that the flow 
effects due to bed-forms penetrated to the free surface of the flow. The 
high velocity flow cell and the relatively high turbulence intensities in 
a region above the trough correspond approximately to the position on the 
surface of natural canals with dune beds at which upwelling of flow is 
often observed. This further confirms the conclusions of Section 4.3.4 


indicating that the bed-forms under investigation were in the dune phase, 
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7.2.3 Pressure Distribution 

The pressure distribution curves were found to be bell-shaped with 
skewness induced by assymmetrical positioning of the maxirmm and minimum 
pressures in relation to the bed-form geometry. The results of pressure 
distribution measurements showed that maximum pressures occurred in the re- 
gion -0.75 £X<-0.55. This corresponds to the mean position of the 
stagnation point established from velocity measurements, The minimum 
pressures were found to occur in the region slightly upstream of ths crest, 
-0.15< ¥< 0.0. This is in general agreement with the findings of 
Raudkivi (1963) who measured pressure distribution over a stabilized 
ripple, but contradicts the results of Vanoni and Hwang (1967) who also 
measured pressures over stabilized bed-forms and found that the minimum 
pressure occurred at or slightly downstream of the crest, 

Downstream from the peak pressure position, the distribution curves 
had inflection points and became concave upward in the region 
~-0,60< X¥ <.0,40, Histograms showing the distribution of locations of the 
maximum and minimum pressures, the maximum pressure gradient, and the max- 
imum bed=-form surface slope are presented in Figure Tele it shows that max- 
imum particle entrainment occurs at the location of the maximum pressure 


gradient, 


7.3 Force Requirements for Particle Detachment 
| 7.301 Statics of Surface Particles 

The mechanics of sediment grain detachment from the surface of a 
bed of similar grains may be investigated by considering the static 
equilibrium of either a discrete grain or 4n assembly of particles in 4 


surface layer, Both approaches were found to be useful in the present 
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investigation. The discrete particle approach provided insight into the 
basic nature of the incipient motion problem For ruteatl work, however, 
the continuum approach was found to be more convenient. 

| In the following paragraphs both approaches are reviewed, compared, 
and an equation for determination of forces required to entrain particles 
from the back faces of active bed-forms is derived. 

In a recent analysis of statics of discrete sediment particles, 
Yalin (1972) showed that the magnitude of an external force required to 
dislodge a stationary particle depends not only upon the weight of the 
particle Wo» but also upon its shape, the geometry of its environment and 
the condition of its surface. 

A typical arrangement of equal diameter spherical grains in a 
two-dimensional situation is illustrated in Figure 7.2(a). Particle "A" 
rests in the hopper formed by two abutting particles "5" and a) Cr. ae ne 
centres of the supporting particles "B" and "C" are assumed to lie on a 
common line parallel to the local surface of the bed-form which makes an 
angle B with the horizontal. 

If the angle between the mean local surface slope of the bed-form 
and the line tangent to particles "A" and "C" at their contact point "1" 
4s V', the normal to the contact point "nl"! makes an angle of V+6 to the 
vertical, For spherical particles of equal size y= 30°. 

The relation between forces at the incipient motion stage for 
particle "A" is shown by the force-polygons in Figures 7.2(b) and (c). 

The body force due to gravity is represented by W,,. which can be estimated 


p 
from equation (7.1): 


Wp = (4/3) r-e(Ps - P) (7.1) 
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where r represents the radius of the particles, g is the acceleration dus 
to gravity, Ps is the density of the particle, and ep is the density of the 
fluid. Ry, is the reaction at contact point "1" and YP represents the 
friction angle between the grains, | 

For a constant Wy the vector of a minimum force pim required to 
disturb the particle equilibrium is shown in Figure 7.2(b), while pFs, in 
Figure 7.<(c), gives the minimum force vector parallel to the slope of the 
bed-forn. 

Clearly, the detachment of the particle may be caused by 4 wide 
variety of mégnitude-.direction combinations of the resultant of the 
peternal impslling forces, 


a few _ tan aS a eat. mise, BY haw = eae oxen : = z se 
The applicatioy of the sine law to the two minimum force conditions 


yields 
pin a Wy sinf a+ 8 ) 7525) 
p's = Wp tana cosp + Wp Sin 6 73) 


where a= W+wW'. 

The equivalent body force diagram and force polygon at threshold 
conditions for a layer of particles resting on a surface inclined 
degrees to the horizontal is shown in Figure 7.3(a) and (b). W is the 
body force duo to gravity of all particles contained in a volume one 
particle diancter thick and one unit in area, and Ww" represents the 
friction angls between the adjacent layers of particles. 


The body force W is given by equation (7.4): 
Woe VC g (Ps -?) = DCY(s -1) (7.4) 


where C is tho fractional volume concentration of particles equal to the 
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ratio of the volume cf solids to the total volume Vir Dis the diameter of 
the particles, g represents the acceleration due to gravity, P, is the 
density of the particles, @ is the density of water, Y is the specific 
weight of water and s is the specific gravity of particles. 

From the force polygon, the external impelling unit force parallel 
to the surface slope at particle threshold, F, is given by equations 


G7e5) or (7.6): 


Sy 
H 


W tan’ cos B +W sin B CRS) 
FPF = DC (s-1)(tan ¥" cosB + sin B) (7.6) 


It is seen from equations (7.3) and (7.5) that the discrete par- 
ticle and surface layer approaches are equivalent provided that ¥" is the 
friction angle which combines the individual particle environmental 
geometry (V ) and surface friction (¥) effects. 

7e3.2 Concentration and Friction 

Equation (7.6) for the minimum force required to entrain the par-~ 
ticles has two terms which indicate thet forces which resist the detach- 
ment arise because of interparticle friction and the direction of dis- 
placement in relation to the gravity field. Both components depend upon 
the thickness D and the concentration of the particles C in the layer 
immediately before entrainment. In addition the friction factor depends 
upon the instantaneous angle of friction between the layers an Both com= 
ponents depend also upon the local surface slope 8 which was the only para- 
meter measured in the expsriments. 

To determine applicable values of C and y" an extensive review of 
literature was carried out, and considerable information related to gran- 


ular particles in bulk volumes was found, However, little has been written 
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about these conditions at the surface of particle beds. Consequently, the 
analysis adopted 4 limit approach; probable maximum and minimum values of 
C and y were estimated from bulk volume considerations and were used to 
determine the limiting bounds for minimum force requirement to entrain the 
particles. 

Theoretical considerations by Graton and Frazer (1935) have shown 
that the volume concentration of equal size spherical particles depends 
upon the packing arrangement of the particles, Maximum and minimum con- 
centrations of C = 0.741 and C = 0.524 are obtained when the spheres are in 
the ideal rhombohedral and cubical arrays respectively. 

| Although in the present experiments spherical particles of nearly 
uniform size distrib 
haphazard as they were deposited by sliding down the lee slope ef the 
bed-forms. 

This type of deposition produces relatively loosely packed granular 
beds, because as the granular solids shear, the mass of particles experige 
ences an increase in bulk volume without change in the volume of solids. 
This well-known phenomenon, termed dilitancy, has been observed in the field 
by Imbrie and Buchanan (1965). They confirmed that avalanched deposits, 
such @S oceur on the lee face of fluvial dunes were loosely packed in re~ 
lation to sheet deposits. 

These observations indicated that the expected concentration of 
particles in the bed-forms of the present investigation should be in the 
lower range of the theoretically possible concentration values. 

é The limits of volume concentraticn of haphazardly packed equal 
spheres were investigated by Scott (1960), In a theoretical study supported 


by experimental observations he determined that the maximum possible con- 
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centration was 0.637 and the minimum possible concentration was 0,601. 

Later, in an experimental study of the effects of various materiel 
properties upon the packing of spherical particles, McRae and Gray (1961) 
found.that 3.08 millimetre diameter spherical particles yielded volume 
concentration of 0,602 when deposited by cascading in air. 

These results were found to correspond closely with the results of 
the average volume concentration of 0.598 of the polystyrene particles 
used in this study as determined from packing of particles for the perneabi- 
lity tests described previousiy in Section 3.3.4. 

From the above observations it was concluded that a value of 0.601 
for the volume concentration of particles in the bed-form could be reason- 
ably assumed. 

For the present experiments, in the absence of better information, 
4% was assumed that the angle of friction between layers of particles on 
the surface of the bed-forms, Y , can be represented by the angle of 
initial yield $4. 

Although in the present problem C and YL" appear as independent 
parameters, past investigations of bulk volume behaviour of granular par- 
ticles have indicated that they are closely interrelated. The relationships 
between C and dj have been extensively reviewed by Allen (1969, 1970a, 
1970b) who concluded, ,from theoretical considerations, that the angle of 
initial yield >; is dependent upon the manner of packing and upon the 
orientation of particles relative to the direction of yield. 

. As numerous ways of particle packing are possible which yield 
different concentrations, Allen (1970a) proposed a semi-theoretical re- 


lationship between C and >: 
tan o>; = 3.596 - 1.884 + A RHP) 














sive 
* ‘ 
7 7 
f 

« » «& ‘ - mag rete Ct le ¢ (sr het Pay "of? Baa 

i ' = Sn ehods lis ne 

o ® ‘ 7 RA'sSO CVSS whew Pit! Tea 

€ 

r . » = 
J 2 < ¢ 7 -* ne “ os 4 ate bo (iiseq or 


- * “ * - @ te - ba 
> 7 f ro; 


ihe erctctee SI * 


———sP* 

y untlosg wort Seateretob aa shade apne 
ye yp 
Ls ratio. ‘ votvotd & eitrbedh "ete a“ 


pvt 4 7" 7 4 ~— rhe ; » a - ff ot x ’ ‘1 oud ~4 a ~ ote ‘ en ¥ wax 
as Te 
, ‘ » 7 <2 = ‘ bh aale ae Dae 
‘ . % " a* - -» oP me > tego » a] 5 SavIanemn * oh cr* : % ae i. 
ra i AGA Shs! os cotton any “T) 


’ ees 


rant! mits nl ,etnentvexe Soesnewd atk —_ 
7 : 


n o nm i ae tenet > . rate 
. J te. ant “Ss + tell nea 


nif 


h Ge ' engety od daa.” Y ,autralcked etd xo Tae 


a 

- . ~ . 7" stun Ve 'é 

, . aa, f w 3 woldory’ Saeaera ‘wld 

= 4 
* 
> . i. im a Seed ™ 
Q ~ "2 - iy A. 2G PJaUvR 
. val 


' aa _ ’ 
aw Pale de + . m ow F ar lf anes wae (+ *tarl4 P ad 
aq Sengolvsiot Of !- .SOT@LeTrIorTal VsiO80 70 * Poa aon eee 


— 
ePOTCE , GOCE) Makin’ y 4 tertvor eed ail niet 

® = : 7 o " >= ms = 
~ Se ofeae edt fade (merodterohs ros {a tes a0) oe ybeks is. ne 


ye 


+ oy pe hie a ated ‘ron ot, He inet toys! 
st : lai ——_ ; 





175 


where A is a dimensionless coefficient to account for real effects such 

as Praacayirieitdon between particles, This relationship implies correct- 
ly that when A = 0 and C is the minimum possible concentration 0. 524, by 
equals zsro. For haphazardly packed equal spheres with A = 0, application 
-of equation (7.7) with maximum and minimum possible concentrations of 
0.637 and 0.601 determined by Scott (1960) yields ob; = 22.1° and 

o; = 155° respectively. Experimental results, in which 0.59 to 0.84 
millimetre size range glass heads were employed (Allen, 1970a), showed 
that A was 0.220 only when C was greater than approximately 0.610, For 

¢ one less than 0.610, the experimental data indicated a very rapid de- 
ae in the d; value with reduction in C, and at C = 0,601 O3 was 
considerably less than 26° predicted by equation (7.7). Ltts: obvious that 
this relationship was not applicable in the present situation, where C was 
established to be 0.601. | 

The value of $4, which represents the angle ¥, was therefore esti- 
mated from experimental considerations. 

The tilting tests described in Section 3.3.5 showed that depending 
upon the mode of initial deposition the angle of initial yield >, was 
either 29° or 31°. The smaller angle of 29° was obtained when avalanched 
deposits of particles were tilted in the opposite direction. And since 
this condition of deposition approximated the conditions of particles in a 
bed-form, >; = 29° was assumed to represent Y" in the subsequent calculations. 

The C and ¥" values of 0.601 and 29° respectively were assumed to 
represent the upper bound of conditions on the surface of the bed-forms. 

For the lower limit, the following observations were considered. 

The concentration C of particles in the surface layer may be substantially 


lower than the particle concentration in the bed-forms because surface 
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layers are subject to dilatation effects arising as a result of shearing ~ 
action by the flow and the impact by other particles, Visual observations. 
such as particle oscillation prior to entrainment and the loose appearance 
of the bed surface confirm the tendency toward a lower C value. Further. 
more, Bagnold (1955 ) has indicated that to allow grain movement to teke 
place at all C must be approximately 0.5, which is the lower limit of the 
theoretical minimum concentration of spherical particles. 

As discussed previously, the concentration of particles is related 
to the angle of initial yield Die To achieve a lower volume concentration 
of particles in ths surface layer while maintaining a constant areal con- 
centration of equal number of particles per unit area, the effective thick- 
ness of the surface layer must be increased by increasing the distance of 
separation between the two upper layers of particles. It would eocean 
therefore that a finite amount of energy must be expended to produce di- 
lation and a slight uplift or rearrangement of surface particles into a 
cubical array. This condition corresponds to an unstable particle arrange- 
ment and the angle of initial yield Os becomes zero. Therefore for the 
lower limit, the concentration of particles in the surface layer was 
assumed to be 0.524 equal to the theoretical minimum for equal spheres 
with Y= 0. 

To compute minimum force requirements to detach particles from 
various locations along the surface of the ped-forms using equation (7.5), 
Y was assumed to be 62.4) pounds per cubic foot, and the specific gravity of 
particles s was adjusted in each test for temperature effects, The local 
value of the angle was determined from the bed-form geometry measurements 


assuming linear profile between limits of individual segments 0.1 L long each. 
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In all calculations D = Deo of 1.4 millimetres was used. This may 
be questioned on the grounds that sliding on the lee face of fae bed-forns 
| produced sorting of the particles in the bed-form, As sorchs by Bagnold 
(1954), coarser particles accumulate near the bottom of the bed-forms. 
Hence, upon re-exposure on the upstream face of the bed-form the mean 
diameter of particles Dso in the zone upstream of X = ~0.50 was likely to 
be larger than in the region -0.50 < X < 0.0. In accordance with equation 
(7.6), the minimum unit force required to entrain the particles from the 
upper half of the bed-forms should have been less then from the lowsr half.. 
However, aS in these experiments the particle gradation was very uniform, 
the errors introduced by neglecting the sorting effect upon local D ee 
would have been of relatively small magnitude in comparison with errors 


in estimating such factors as concentration or friction. 


7.4.1 Boundary Shear Stresses 

No reliable methods for aeteraining boundary shear stresses on the 
surface of active bed-forms have been established to-date, The problem is 
exceedingly difficult because the combined effects of pressure gradients, 
movement of solid particles, curved boundaries, and inflow and outflow 
from the porous bed must be considered simultaneously in establishing the 
proper relationships. In the absence of a better approach, the shear stress 
distributions on the upstream face of the bed-forms in the present investi- 
gation were estimated using plane rigid boundary techniques. 


The application of rigid boundary techniques to determine skin 


resistance to flow over alluvial bedeforme has been attempted by Einstein 


and Barbarossa (1952), Simons and Richardson (1966), Vanoni and Hwang (1967), 


and others, Generally these investigators assumed that the local shear 
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stresses can be-estimated-from some type of an universal logarithmic velo. 
city relation. 

The mean dimensionless velocity profiles Te presented in Figure 
6.13 disclosed that between X = -0,50 and X = -0,10 and near the Ee ee 
boundary the velocities were logarithmically distributed. Therefore, to 
determine skin friction distribution characteristics the law of the wall 
relationships was employed. Preliminsry considerations showed that at 
y/d = 0.15 the shear velocity (ux) Reynold's cmeersenay/e in the region 
downstream from the reattachment point were in the order of 100. The 
thickness of the viscous sublayer §' calculated from 11.6.¥/u,, was found 
me bs approximately 0.003 feet. Hence the ratio Deo/s' of the character- 
istic particle size to the thickness o 
to be about 1.5. This represents conditions of flow over a hydrodynemically 
smooth surface in which the surface roughness elements are just beginning 
to have effect. Therefore, for surface shear stress determination the 


smooth boundary law of the wall relation given by equation (7,8) was 


employed: 
u i yu 
ES eee os: (7.8) 
Ux K Vv 


where u is the mean velocity of flow at depth y, v is the kinematic vis- 
cosity, B is an experimentally determined constant, and x represents the 
yon Karman's constant. According to Goma and Gelhar (1968), the value of 
k = 0,40 was consistent with observations of flow in pipes with porous 
walls, It was therefore assumed to be applicable also in the present 
experiments. 


The relationship to determine boundary shear stresses was obtained 


by rearranging equation 7.8 into 
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us (5.75u,)log y +C (7.9) 


whers c'4s a constant tern not containing either u or y. Therefore 
(5.75u, ) represents the slope of the velocity profile in au v.s. loz y 
plot. - The method is constrained, however, to logarithmic velocity distri- 
butions at depths T < 0.15. 

The information required to establish velocity distribution slope 
for each shear stress evaluation was obtained from experimentally determined 
data, The mean dimensionless velocity profiles 7 VeSe x » together with - 
mean depth information given by the individual bed-form geometry data and 
the measured flow discharge in the flume were used to reconstruct the local 
u vs. log y profile. From this the slope of the velocity profile 
du/d(log y) between y/d = 0.0625 and 0.125 could be determined. The Local 


shear was then computed from equation (7.10) 


% = 4 are aCe & - (7-10) 
The results of shear stress evaluation in all segments downstream 
from X = ~0.50 above all selected bed-forms are presented in Figure 7 aka 
From physical considerations, and as indicated by the results of Laursen 
et al, (1962) and Raudkivt (1963), the shear stress at the reattachment 
position, considered fron velocity profile data to be at X = -0.65, was 
Meeened to be zero. At X = -0.55 shear stresses were not calculated 
because the measured velocity profile as shown in Figure 6.13 didnot 
appear to be logarithmic. The reason for this anomaly at X = -0.55 is 
not clear. However, this location corresponds generally with the location 


of the maximum pressure gradient, the maximum bed-form surface slope, and 


the local suppression of turbulence 45 discussed in Section 7.2.2. 
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Figure 7.44. Boundary Shear Distribution. 
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| The results of the average shear stress distributions on the 
boundary of the bed-forms show that the shear stresses cee ene iiee 
from the reattachment position downstream to X = -0.45, and remain constant 
to X = -0.15. However, downstream from X= ~0,15 the shear stresses dimin- 
ish rapidly as indicated by the relatively low value of shear at X = -0.05 
and appear to be zero at the crest. These results are clearly not in agres- 
ment with the boundary shear stress fg ates ss ee established on stabilized 
bed-forms by Laursen et al. (1962) and Raudkivi (1963) who found that the 
shear stresses increased continuously from the reattachment point downe= 
stream to the crest of the bed-forms. 

The discrepancies between the ieee findings neers results 
obtained from measurements on rigid boundary bed-forms were attributed to 
the effects of particle movement along the bed and the possible outflow 
from the bed immediately upstream from a crest. _ 

In the present experiments boundary shear stresses were established 
from velocity profile data which assumed boundary of the bed at the upper 
surface of the sediment particles. Because of bed particle movement the 
assumed boundary, as indicated in Figure 6.12, had a finite average veloe 
city, and furthermore was located above the stationary particles. It 
appears that the velocity of the fluid in the layer adjacent to the parti- 
cles was increased because of particle movement. This reduced the velocity 
gradient and, therefore, the estimated shear. 

On rigid boundary bed-forns, the shear was found to increase because 
as the boundary layer developed in the downstream direction from the re- 
attachment point, the velocity gradients near the boundary increased. 
Boundary layer developed on 4 granular boundary also. However, it is possible 


that the movement of the entrained particles, increasing in velocity in the 
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downstream direction, increased sufficiently the local velocity of fluid 
above the particles and effectively counteracted the increasing velocity 
gradients due to boundary layer development. 

Experimental evidence suggested that small amounts of outflow 
from the bed into the flow occurred in the region upstream of the crest. 
Although pressure measurements did not disclose any vertical gradients, 
horizontal pressure distribution measurements showed that minimum pressures 
occurred immedistely upstream of the crest. Therefore, pressure gradients 
would have produced seepage toward the low pressure zone from both sides, 
Evidence to support this observation was provided by the direction of 
particle movement near the crest on the lee side of the bed-form imne- 
diately prior to deposition, as diseussed in Section 6.5.1, and also by 
the extremely steep lee-side slopes of the bed-forms as discussed in 
Section 7.2.1. As the flume bottom was impermeable,to maintain continuity, 
the outflow must have occurred through the interface of the bed-form up- 
stream of the crest. 

A summary of qualitative effects of fluid injection into the flow 
upon channel boundary leyer has been presented by Watters and Rao (1971). 
According to this summary injection would have reduced the friction drag 
upon the boundary as was confirmed by the reduction of shear estimated 
from the experimental evidence. 

Although the shear measurement results are subject to the limi- 
tations imposed by the measurement and evaluation technique, they never= 
theless point out some possible misleading interpretation of results 
obtained from measurements on simuleted bed-forms. Especially significant 
is the present observation that shear does not increase but decreases 


near the crests of the bed-forms. This seems to be reasonable because no 
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or very little particle entrainment near the crest takes place. If the 
shear were to increase continuously toward the crest, as indicated by 
Raudkivi (1963), more particles would have been entrained from the crest 
position, and the bed-form would have not retained its shape. 

7e4+e2 Forces Due to Pressure Gradients 

The results of pressure measurements, presented in Figures 6,21 
(a) to (h), have shown that pressure gradients were established through- 
out the mass of granular particles in the bed-forms, The characteristics 
of these pressure distributions were found to agree qualitatively with 
pressure distributions established on the surfaces of stabilized bed- 
forms as reported by Laursen et al. (1962) and Reudkied (1963). On this 
basis it was assumed that the measured pressure characteristics were re- 
presentative of conditions in the layer of stationary particles on the 
surface of the bedforms in the present expsriments. 

When pressure gradient exists in the fluid of the pores of granu- 
lar materials, a dynamic force in the direction of decreasing pressure, 
4s exerted upon the stationary particles (Taylor, 1948). In the layer of 
stationary particles resting on the surface of the bed-forms the distri- 
bution of unit force, or stress, Up per unit area due to pressure gradients 


was estimated by equation (7.11): 
tp = hy ey, Deo cos B (7.11) 


where h, represents the local hydraulic pressure gredient in inches of 

water per inch length of bed-form, 7 is the specific weight of water, Deg 
4s the characteristic particle diameter, and 6 represents the local slope 
of the bed-form surface, The pressure gradient was determined from the 


mean pressure distribution characteristics by estimating the slope of the 
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pressure distribution curve at the mid-point of each 0.1 L segment. The 
‘surface slope angle B was determined assuming linear profile between the 
limits of each segment. For all calculations Deg was assumed to be 1.4 
millimetres. 

The distribution of the pressure shear stress Uy, upon the particles 
in a surface layer of the bed-forms downstream from the reattachment point 
are presented in Figure 7.5. 

In the majority of bed-forms investigated, the maxcinmum surface 
slope angle 6 was less than 10°. The reduction in the local value of 155 


due to bed surface inclination was therefore less than two per cent, 


7.5 Relative Force Contribution 
7e5el Results 
_ The relative contributions of the average boundary shear stresses 
due to fluid flowly and pressure gradients Tp on the particles in the 
surface layers of the bed-forms downstream from the reattachnent point were 
determined for the two extreme conditions of particle concentration and 
friction assumed in Section 7.302. 

The percentage contribution of each of the two stresses to particle 
entrainment and their distribution on each bed-form were computed by di- 
viding the local values of Uy or Tp in each segment by the total force 
requirement F determined by equation 7.6. 

The results for the upper limit condition with C = 0.601 and 
py 29° are presented in Figure 7,6 and show that the force contribution 
to the entrainment of particles by flow shear increases rapidly in the 
downstream direction from zero at the reattachment point to approximately 


X= -0.45, gradually thereafter to 4 maximum at X = -0.25 and decreases 
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to zero at the crest. The maximum relative unit force contribution due 
to flow shear occurred at X = -0.25 and was found to range from 27 to 78 
per cent depending mainly upon the mean flow velocity. 

The maximum unit force contribution by the pressure gradients 
Peerred in the region between X = ~0.60 and X = ~0.45, and was found to 
be approximately six per cent in nine out of the thirteen bed-forms, In 
the rem#ining four bed-forms this value ranged from four to eight per cent. 

Obviously, the maximum pressure gradient en eccurred in 
the region of increasing flow shear, This tended to equalize the dis- 
tribution of the combined flow shear and pressure gradient stresses and 
in @ number of bed-forms the combined contribution between X = -0,50 and 
X = -0.20 was found to be relatively constant. 

As shown in Figure 7.6, the sum of the mean flow shear and preS- 
sure predient stress contributions on the upstream face of the bed-forms 
between X = -0.50 and X = -0.15 provided from 30 to 80 per cent of the 
total required stress to entrain the particles. 

Under the lower limit conditions of C = 0.524 and Ue 0, which 
in effect assumed no friction between the particles, the total resisting 
force throughout most of the bed-form face downstream from the reattach- 
ment point was found to be less than the estimated flow shear T,. There- 
fore, only the relative unit force contribution by the pressure gradients 
were estimated and summarized in Figure Tele 

The results show that the maximum relative contribution by the 
mean pressure gradient occurred also in the region between X = -0.60 and 
o -0.45, and in different bed-forms ranged from 30 to 80 per cent. This 
range corresponds in magnitude with the sum total of the relative contri- 


butions by flow shear and pressure gradient stresses under the assumed 
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Figure 7.7. Relative Contribution of Pressure Gradient Shear Stress 
with Frictionless Surface Layer Conditions. 
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-ypper limit conditions, However, possibly because of considerable data 
scatter, no direct correspondence between the sum total of flow shear 
and pressure gradient stress contributions under the Biv concentra= 
tion and friction and the pressure gradient stress under the frictionless 
conditions on individual bed-forms could be conclusively established. 

The maximum relative pressure gradient contribution under fric- 
tionless conditions.and the corresponding maximum estimated local rate 
of particle entrainnent for twelve bedforms were compared in Figure 7086 
The comparison shows thet on eight bed-forms the rate of particle entrain- 
ment increased with increased relative contribution due to pressure graq 
dient stresses. Data points from the remaining four bed-forms did not 
follow the same trend and from the availsble information it was not 
possible to explain the reasons for the discrepancy. 

75.2 Discussion 

Although every attempt was made to obtain accurate experimental 
data, ths results presented in the previous section are generally of 
qualitative nature, A number of assumptions were employed in computing 
the relative contributions of the mean flow shear and the mean pressure 
gradient stresses for the detachment of sediment particles from the sur- 
face of active bed-forms. Rigid boundary techniques were used for the 
determination of boundary shear stresses from mean velocity profiles 
which themselves were developed on the assumption of similarity at 
corresponding positions relative to the geometry of the bed-forms. 
Secondly, with the exception of the local slope of the bed-form surface, 
actual conditions in the surface leyer of particles, necessary for the 


determination of forces resisting particle entrainment, were not experi- 
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mentally obtainable. Hence, the mean measured stresses could not be 
compared with the measured resisting forces, but were compared with the 
calculated unit forces based upon two extreme surface layer conditions 
immediately prior to particle detachment: (1) maximum concentration and 
friction corresponding to assumed conditions in the bulk of particles 
deep within the bed-form, and (2) minimum concentration and no friction 
between particles corresponding to theoretically unstable arrangement 

- of equal spheres, The comparisons were made also considering local drag 
forces parallel to the surface of the bed-forms because the nature of 
experimental data provided no information from which lift forces could 
be evaluated, 

Nevertheless, even under the above assumptive restrictions, the 
qualitative results of this investigation show @ number of previously 
not established facts about the entrainment of sediment particles on the 
backs of the bed-forns,. 

The relative contribution results obtained under the maximum con- 
centration and friction conditions show that the combined drag by the 
local mean flow shear and the mean pressure gradient stresses was not 
adequate to entrain the particles from the surface of the bed-forms,. 
Furthermore, the distribution characteristics of the relative unit forces 
showed that entrainment of particles was occurring from positive bed-form 
surface slope locations where the mean flow shear and pressure gradient 
stresses were calculated to be zero. 

These observations suggested, first, that other force systems 
were acting upon the particles on the surface of the bed-forms, Among 
these may be the hydrodynamic lift as considered by Einstein and El-Samni 


(1949) and the effects of turbulence fluctuations. However, lift by 


“ om ; 
- i, - a a 7% 
we Y »o fa ga we 00 > 180Ne 













wo hs 2 | asS sxe ow? noqn beead sanrt 


mints (S$) bas wesetabed eft 
odd od yatbaoquertoo eolotirag 
a 
9 


txage> eff .woredqe. 
ee ai | 7 

~ ‘ . er > Ps pe - 

At lo eoetrre odd of Leliorag we 


aed 


; ite 
a 2 Lead nm? aind I6 sttnebeiedll 


Be 


® sit Jgods etoal betatiay 


Drtot-be 0 a 


vi soltudistace avis Ses oat 
oe 
- 
f s< bd 7 = > mar. es. 
ef? Jad da enol? ifneo mole peat pee tia 
ny le _ 


PP 
Be . sgnegq Gaea ed? bas aoe — a 
° T othe At 
| s ef? sot? eelotiveg -add 


a 


i 
Rens jiar evitela: efS to tolstefrevortado aotse pd iy stiabo oa 






ec te trea out: =~ ro sad 
sir 






Ly , ae 





daethass eieteng bas qeds wotl 
pious 


203 


itself cannot entrain the particles because translation of the bed-forms 
requires that the particles be displaced in the direction of flow. There- 
fore, drag forces must be present. The experimental results similarly 
suggested that turbulent velocity fluctuations were not directly respon- 
sible for the entrainment. The maximum turbulence intensity of the hori- 
pontal velocity component was observed to occur in the trough region of 

| the bed-form which is the region of zero particle entrainment. 

In addition to other force systems, the calculated local relative 
contributions by the mean flow shear and mean pressure gradient stresses 
may also have been influenced by the concentration and friction values 
used in the calculations. Equation 7.6 shows that the minimum drag re- 
quired for particle layer entrainment has to overcome the friction between 
the particles and the gravity force component which depends upon the 
local bed-form surface slope. The mass of particles in the surface layer 
4s governed by the concentration. 

In evaluating the unit force requirements to detach the particles, 
the friction angle W and concentration C in the surface layer above a 
bed-form were assumed to be constant. This may not be true as visual 
observations and physical considerations suggested that both values may 
have been affected by the action of the flow, turbulence or particle in- 
pact even before the entrainment. Generally the mass of particles on the 
surface of the bed-forms was loose indicating that the assumed concentra 
tion of 0.601 was too high. In some regions of the bed-form, especially 
in the vicinity of the reattachment point, particles were observed to 
oscillate above their supporting particles suggesting that at least the 
static friction between particles had been overcome, In regions of trans- 


port, the impact of moving particles upon the stationary particles 
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exerted an equivalent shearing action which dilated the stationary mass 
of particles and reduced the concentration, These observations suggested 
that the conditions in the stationary surface layer of particles may be 
substantially altered prior to actual particle entrainment, 

In relation to these observations, it is Significant that quali- 
tatively the pressure gradient distribution characteristics were found 
to correspond more closely with the particle entrainment distribution 
curves than other parameters, The maximum entrainment rates were found 
to occur at the location of maximum pressure gradients, Zero or very 
small entrainment occurred immediately upstream of the crest of the bed- 
forms, where the pressure gradient was zero but where the shearing action 
by particle impact was likely to be most intensive and dilation most 
complete, Near the naa LEStmenceeaints although the mean pressure gra- 
dient was zero, the finite entrainment could be attributed to pressure 
fluctuations, 

These observations suggest that pressures and their gradients 
may be the basic variables governing the entrainment of particles from 
the surface of the bed-forms. This is in agreement with the comments of 
Einstein and El-Samni (1949) regarding statistical description of tur- 
bulence in the vicinity of a wall and Luque's (1967) theory of stability 
of erosive soils, However, if pressure gradients and pressure fluctuations 
actually are the basic variables, then boundary shear and velocity fluc- 
tuations are relegated to other functions--perhaps to the dilation of the 


surface layer particles, 
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CHAPTER VIII 
SUMMARY, CONCLUSIONS, AND RECOMMENDATIONS 


8.1 Sumary 

The present experimental investigation was undertaken to obtain 
a better insight into the mechanics of bed-forms in ailuvial channels, 

A distinctly new experimental approach was adopted for the 
study. The measurements of various flow and bed-form parameters were 
made using actively translating bed-forms with bed-load mode of sediment 
transport. This differed fundamentally from similar studies carried out 
by other investigators in the past who employed simulated, non-translating 
bed-form shapese | 

The development of experimental techniques for use in active 
sediment transport and bed-form translation environment formed a 
significant part of the investigation, All experiments aye made in a 
Laboratory flume with spherical, light weight (polystyrene) particles of 
nearly uniform size gradation. Measurements were made of the pressure 
distribution characteristics within the mass of granular particles in 
the bed-forms, of the velocity profiles in the flow above the bed-forms, 
and of the bed-form ceometry. 

All instrumentation probes were stationary.in the flume and all 
instantaneous observations of the various parameters were related to the 
geometry of the bed-forms, Also, all probes were of sufficiently small 


size and disturbed neither the local flow field nor the sediment 


boundary, 
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Pressure distribution characteristics were establishe 
differential pressure gages, Difficulties were experienced in the 
measurement of pressures because of vsry small pressure variation above 
light weight sediment bed-forms, pressure fluctuations dus to surfece 
wave activity and turbulence, changing roughness in the flume due to 
natural development of bed-forms, and the loss of water from the re- 
circulating system of the flume. Special techniques were developed and 
successfully employed to minimize these effects, 

Flow velocity measurements were made using hydrogen bubble 
techniqus. This technique was found to be satisfactory because the 
inclusion of the small diameter bubble generating wires had absolutely no 
effects upon the flow or the sediment bed, It provided the information 
for establishment of instantaneous velocity profiles throughout the full 
depth of flow, and rescived flow in upstream and downstream directions. 
The major disadvantages of the hydrogen bubble method were the large 
number of instantaneous observations necessary to establish mean velocity 
profiles, the extremely time consuming processes of data reduction, and 
difficulties in data interpretation near the sediment bed. 

Measurement of bed-form profile characteristics were made using 
& series of small diameter wire sceles installed vertically along the 
centreline of the flume. 

Bed-forms wore naturally formed, two-dimensional and in the dune 
phase. Because of secondary circulation effects which tended to distort 
the bed-forms inte irregular features, two-dimensionality of the bed- 
forms was maintained by manually advancing the retarded crests of the 
bed-forms near the flume sidewalls. 


Major difficulties were experienced in obtaining constant bed~ 
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form hoights and lengths during thotr passage through the observation 
area, Statistical techniques were used to select seventeen of the 


fifty-five bed-forms observed in the experiments, 


Ban Conclusions 

The conelusions reached from te results of this investigation 
were as follows: | 

(1) Pressure gradient induced seepage forces contribute a 
finite force upon the stationary bed particles in the surface layer of 
the bedeforiit. 

(2) The maxima pressure on the bed-forms eccurs at the re- 
attachment point, and the minimum at or slightly upstream of the crest. 

(3) The lecation of the maximum mean pressure gradient was 

approximately one tenth of the bed-form wevelength downstream from the 
reattachnent point. 

(%) The location of the maxinum mean pressure gradient corre= 
Ponds with the position of the maxirmm bed-form surface slope &nd 
with tho region of maximum particle entrainnent. 

(5) Tho minimum pressure gradient occurs slightly upstream 
of the crest of the bed-forms and corresponds to the location of minimum 
particle entrainment from the upstream face of the bed-forms. 

(6) The mean boundary shear, estimated from mean velocity pro- 
files by rigid boundary techniques was found to increase downstream from 
the reattachment point to a point approximately one quarter of the wavee 
length upstream of the crest and to decrease to sero at the crest. 

(7) The mean local velocity of flow occurs between 0.3 end 0.4 


of the local depth above the bed-form surface. 
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(8) The maximum horizontal velocity turbulence intensity 
occurs in the trough region and approximately 0.10 above the surface 
ef the bed-forn. 
(9) The horizontal velocity turbulence intensity appears to be 
suppressed in the vicinity of the maximum bed-form surface slope. 
_* (10) From the trough region, turbulence effects penetrate to 
the surface of the flow. 

(11) The relative contributions of the mean flow shear and pres- 
sure gradient stresses upon the particles contribute only a part of the 
total required forces for particle entrainment. 

(12) The mean pressure gradient distribution above bed-form 
surface downstream from the reattachment point corresponds with the 
particle entrainment distribution characteristics except between the 


reattachment and maximum entrainment locations. 


8.3 Recommendations 

Further studies are necessary to establish quantitatively the 
contributions of pressure gradients and pressure fluctuations, boundary 
shear stresses and turbulence upon the entrainment of particles from the 
surface of nade rovns: Specifically the following areas should be 
investigated: 

(1) The determination of boundary shear stresses on porous 
boundaries with pressure gradients, surface curvature and movement of 
discrete particles. 

(2) The effects of flow shear and turbulence upon entrainment 
and/or dilation of surface particles, including particle resonance to 


fluctuations and particle gradation effects. 
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(3) Pressure distribution, pressure fluctuation and seepage 
characteristics within and near the surfece of bedforms, 

It is recommended that, because of difficulties in pressure 
measurement, heavier bedematerials be used in future studies, If 
hydrogen bubble technique is used for velocity measurement, the bed 
materials should also be of dark colour and non-reflecting. Al] instru. 


mentation should be automatized. 
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| TABLE A.1 
SEDIMENT PARTICLE SPECIFIC GRAVITY 











Test Temperature %¢ Specific Gravity 
ng 25.6 1.0427 

2 2526 1.0366 

3 Zl. 5 1.0430 

& 24.5 1, 0425 

5 25,7 1.0383 

6 2507 1.0415 








TABLE A,2 
SEDIMENT PARTICLE SPECIFIC GRAVITY DISTRIBUTION 


ACNE CS SEP EEE LT LSI TEN ELIA IIT IIE LOE BIE IE AAT AES TT 
SRRCREED Sry ares 














Sdicads = 














RAST 





Ranges of Per Cent by Weight 
Specific Gravity Test #1 Test #2 Test #3 Test #+ Average 
1,020 - 1,025 0 0 0 1 0025 
1,025 = 1,030 u 0 2 0 0.75 
1.030 » 1.035 6 iil 8 8 Oseo 
1,035 - 1,04 7 19 ay 23 19.00 
1.040 = 1.045 69 58 57 6h, 61.50 
1.045 ~ 1.050 18 9 at 0 7.00 
1,050 =~ 1.055 a 3 5 3 3209 
1.055 ~ 1.060 0 0 0 1 0.25 


eqecerner este zu I SH OE EA RE IO SES TTT IESE TED CE TE SE ES a TT 


TABLE A.3 
SUMMARY OF GRAIN SIZE DISTRIBUTION TESTS 





Size of Opening Per Cent Finer by Weight 


Millimotros Test #1 Test #2 Test #3 
nee 
4,760 100.00 100,00 100,00 
2.000 98.45 98,94 98. 50 
0.840 5557 0.87 0.92 
0.591 0.25 0.39 0.23 
0.420 0,06 0,07 0.10 
0.149 0,00 0,00 0.90 
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TABLE A.4 
PERMEABILITY AND VOLUME CONCENTRATION DATA 





GRISLY SALTER LR PTT NTS IEE HD 


Coefficient of 




















Volume 








Sample Permeability* Porosity Concentration 
c cm/sec 
A 1.16 x 107! 0.377 0,620 
B 1.09 x 107! 0.377 0,620 
C 1.19 x 1074 0.459 0. 540 
D 1.32 x 1074 0.399 0.60% 
E 1.30 x 1o~1 0.399 0,604 
Rvarace Loa a yo~l 0,402 0.598 





* At 20°C, 
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PROCEDURE A,1 
METHOD TO DETERMINE DISTRIBUTION OF SPECIFIC 
GRAVITIES OF SEDIMENT MATERIAL 

(1) Place 50 grams of dry particles into distilled water for 24 
hours to climinate air bubbles adhering to the surface of the particles. 

(2) In risers used for soiis hydromster analysis prepare eight 
one litre solutions of differont specific gravities by dissolving different 
amounts of sodium phosphete detergent (Calgon) in water, 

(3) Remove the sample from the distilled water and place in the 
sodium phosphate solution with the highest specific gravity. Stir the 
nixture vigorously for one minute, and remove particles which rise to the 
surface with a special mesh spoon, Allow the particles to drain and place 
them in a cylinder containing solution with the highest specific gravity. 
Repeat the process to the last solution with the lowest specific gravity. 

(4) Using a hydrometer determine the spscific gravity of the 
solutions with the remeining sediment particles in each cylindrical 
flask, pee sayecas! + moasure the temperature of the solution, 

(5) Remove the sediment particles remaining in each sedimentation 


cylinder, rinse with distilled water, dry at 30°C temperature, and woigh. 
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LIST OF BQUIPHENT FOR HYDROGEN BUBBLE PRODUCTION, 
CONTROL, AND PHOTOGRAPHY 
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1. Electrical Ecauipnent 
(1) Power Supply: 
Regulated Power Supply Unit, Lambda Corp., Model C-1581M, 
(2) Control Equipment: | 
(a) Pulse generator; Hewlet-Packard, Model 214A 
(b) Frequency oscillator: Krohn-Hite Corp., Medel 400-C 
(c) Relay, mercury wetted type; Potter and Brumfield Coes 
Model JM 311531 


2. Photographic Equipment 
(1) Camera: 
Nikon F with f 3.5/55 mm macroelens,. 
(2) Film: 
Kodax 2475 Recording Typs rated at 1000 ASA, 
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SUMMARY OF STATISTICAL BED.FORM STABILITY 
TEST RESULTS 
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